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Strange Quark Mass Determination from 
Sum Rules for Hadronic r-Decays 

Abstract 

We discuss the ratio of hadronic to leptonic r-decays, that can be expanded in 
an operator product expansion. The sensitivity to the strange mass is increased, 
if only the flavor-breaking difference of strange to non-strange currents is consid- 
ered, which leaves us with terms that are at leading order proportional to m 2 s . 
In this sum rule the scalar channel exhibits a very bad convergence behavior 
in the perturbation series, that we control by replacing the scalar OPE with a 
phenomenological ansatz for the scalar and pseudoscalar spectral functions. The 
spectral ansatz is compared to the theoretical (OPE) expressions and we show 
that the results agree within the uncertainties and that the uncertainties of our 
ansatz are smaller than those of the theoretical expressions. This allows us to con- 
siderably reduce the uncertainties in comparison to previous m s determinations 
from r-decays. Our final results are m s (2 GeV) = 117±17 MeV for the unitarity 
fit of \V US \ and m s (2 GeV) = 103±17 MeV for the Leutwyler-Roos value of \V US \. 



Strange-Quark Massenbestimmung aus 
Summenregeln fur Hadronische r-Zerfalle 

Zusammenfassung 

Wir untersuchen das Verhaltnis von hadronischen zu leptonischen r-Zerfallen, 
das in einer Operator Produkt Entwicklung entwickelt werden kann. Um die 
Sensititivitat auf die Masse des Strange-Quark zu erhohen, wird nur die flavor 
brechende Differenz aus ud and us Stromen betrachtet, die in fiihrender Ordnung 
proportional zu m 2 ist. Der skalare Kanal der so erhaltenen Summenregel zeigt 
ein extrem schlechtes Konvergenzverhalten in der Storungsreihe, das wir umge- 
hen, indem wir die skalare OPE durch einen phenomenologischen Ansatz fur 
die skalaren und pseudoskalaren Spektralfunktion ersetzen. Unser Ansatz wird 
daraufhin mit den theoretischen (OPE) Ausdriicken verglichen, und wir zeigen, 
daB die Ergebnisse im Rahmen der Unsicherheiten ubereinstimmen und die Un- 
sicherheiten unseres Ansatzes kleiner sind als die der theoretischen Seite, was uns 
erlaubt, die Fehler im Vergleich zu fruheren m s Bestimmungen aus r-Zerfallen 
signifikant zu reduzieren. Unsere Ergebnisse sind m s (2 GeV) = 117±17 MeV fur 
den Unitaritatsfit von |T4s| und m s {2 GeV) = 103±17 MeV fur den Leutwyler- 
Roos Wert von \V US \ . 
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Chapter 1 
Introduction 



We have a simple, definite theory that is supposed to 
explain all the properties of protons and neutrons, yet 
we can't calculate anything with it because the mathe- 
matics is too hard for us. 
(Richard P. Feynman, QED) 

If one were to ask a modern day elementary particle physicist about his assess- 
ment on the status of the field, the answer would almost certainly be ambivalent: 
on the one hand we have an extremely well tested theory, the standard model 
(SM), that seems to accurately describe all observed phenomena. On the other 
hand this resistance against falsification has, to an extent, turned into a problem: 
there are quite a few "unsatisfactory" elements in the theory, that demand further 
investigation or should not be present in a fundamental theory: the most glaringly 
obvious is of course the fact that gravitation is not included at all; furthermore 
we do not understand the Higgs sector, and one would hope for a smaller number 
of fundamental parameters that have to be fixed by experiment. In the simplest 
version of the standard model there are 18 free parameters, among these param- 
eters the quark masses are the ones less precisely known, along with the Higgs 
mass and the CKM matrix parameters. These uncertainties in the quark masses 
are mainly due to the fact that quarks (except for the top, that shall be excluded 
from the following discussion) are not observed in nature as free particles, and 
their determination is then linked to another, more technical, problem of the SM: 
the theory that is supposed to describe all effects of quarks and their "binding 
forces", the gluons, Quantum Chromodynamics (QCD) pQ, while perfectly tested 
and calculable at high energies due to asymptotic freedom (the smallness of the 
coupling constant) fails at low energies when the QCD coupling becomes large 
and perturbation theory breaks down. This is known as confinement and leads to 
the fact that the observed degrees of freedom are the hadrons not the fundamen- 
tal quarks and gluons. It is then clear that, in order to determine quark masses, 
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nonperturbative methods will have to be used. 
Today, the most important of these are: 

• Chiral Perturbation Theory: Introduced by Weinberg |2j, Gasser and 
Leutwyler [311], Chiral Perturbation Theory is an effective field theory that 
begins from the notion that the QCD Lagrangian is symmetric under chiral 
transformation in the limit of massless quarks. The lowest lying meson 
octet is seen as the Goldstone bosons of the spontaneous breaking of this 
chiral symmetry. Additionally, the symmetry is also broken explicitely by 
the quark masses and the idea of the theory is to calculate for example 
mesonic scattering by constructing a perturbation series in small momenta 
and quark masses. Furthermore, (and for our purposes more interestingly) 
one can also calculate mass ratios j^]. An absolute scale, however, can only 
be determined by another method. 

• QCD Sum Rules: Introduced by Shifman, Vainshtein and Zakharov 
[3 Ej(SVZ), sum rules are up to now the most successful method for the 
determination of quark masses. The basic idea is to relate hadronic quan- 
tities to fundamental QCD parameters by using dispersion relations and 
the operator product expansion and to include nonperturbative effects by 
defining operator vacuum expectation values as condensates. It is inter- 
esting to note that originally sum rules were used to determine hadronic 
properties by giving the fundamental parameters as an input; this has been 
turned around when experimental precision on the hadronic side improved 
over the years. 

• Lattice QCD: Based on the idea of the Wilson Loop [H| it is possible to 
discuss QCD on a lattice instead of in the continuum and this will almost 
certainly lead to the most precise determinations in the future. Up to now 
the calculations have been performed mainly in the so called "quenched" 
approximation that neglects internal fermion loops. Mathematically this 
corresponds to setting the fermion determinant to 1, which leads to a sig- 
nificant improvement in computational speed. It is, however, not entirely 
clear how good this approximation is. 

Further nonperturbative methods that have less impact on quark mass determi- 
nation are phenomenological quark models, in which the quarks are assumed to 
move in an effective potential (examples for this would be the MIT Bag or the 
Isgur Karl model) and expansions in 1/N C , where the number of colors N c is 
assumed to be large. 

In this work we will concentrate on the strange quark mass and determine it 
from QCD sum rules for r-decays. There has recently been some activity in 
both the sum rule CHI HU 021 03] and the lattice [111 [Tj3 EE] community for 
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a determination of this parameter, but the relative uncertainties remain higher 
than for the heavier quarks, while the best determinations for m u and use 
chiral perturbation theory and m s , which causes their uncertainties to depend on 
those of the strange mass. It is noteworthy that m s is not only interesting as a 
fundamental parameter but also shows up in the standard model calculations of 
the directly CP violating parameter e'/e. 

The first chapter of this work deals with the basic theoretical background: we 
begin by very briefly sketching renormalization from the QCD point of view in 
order to motivate the renormalization group equations and the running of cou- 
plings and masses. Additionally we introduce the main ingredients for QCD sum 
rules, the dispersion relations and operator product expansion mentioned above 
as well as the types of sum rules most commonly used today; among these are 
the finite energy sum rules that the rest of this work will be dealing with. 

In the next chapter we show how finite energy sum rules can be applied to 
r-decays in order to determine m s . The advantage of these sum rules is the 
availability of rather precise data from r-decay studies for example at ALEPH 
[T7| ITS! \T§\ I2U] . The OPE for both scalar and vector channel of the sum rule 
are then given, before we go on to examine the convergence of the expressions 
in both channels. We observe an extremely bad convergence in the perturbative 
contribution to the scalar channel, that causes large uncertainties in a strange 
mass determination from the r-sum rule. 

In the third chapter we present the phenomenological ansatz that we make 
to avoid this problem of convergence: we use spectral functions for scalar and 
pseudoscalar us and ud currents. These spectral functions can be considered as 
resonance dominated in both pseudoscalar channels, while the scalar us channel 
has to be treated dynamically. Previously calculated scalar form factors enable 
us to provide such an ansatz. The sum of all phenomenological expressions is 
dominated by the kaon pole and our phenomenological ansatz is thus afflicted 
with smaller uncertainties than the corresponding theoretical expressions. In this 
chapter we also investigate briefly the possibility of using our spectral function 
for an m s determination from pure scalar sum rules. 

Finally, in the last chapters we give the experimental data and the numerical 
results that we find from them. Here we discuss also the use of the input param- 
eter \V US \ and the significant impact on our determination as well the results of 
other determinations. 
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Chapter 2 

Theoretical Background 



There is nothing so practical as a good theory. 
(Richard P. Feynman, QED) 



2.1 The three R's: Regularization, Renormal- 
ization and the Renormalization Group 

2.1.1 Regularization 

A thorough discussion of the complete renormalization program can be found in 
any good textbook on quantum field theory of which there are many, so only 
the most important issues shall be reviewed briefly, emphasizing those methods 
that are relevant to QCD. The general approach for any field theory is to start 
from a Lagrangian density that consists of kinetic and interaction terms. From 
this Lagrangian one can derive Feynman rules that allow transition amplitudes 
to be written as a sum of Feynman graphs. In an ideal theory (which QCD 
is not, we will see why later in this chapter) the coupling constant a is small 
enough to admit a perturbative expansion in a. All is well as long as one is 
satisfied with tree-level diagrams that do not involve loops. It is, however, in the 
human nature not to be content with this approximation and to strive for further 
precision. Unfortunately, a naive computation of a diagram like in Fig. 12.11 that 
contains loops turns out to be divergent, because it involves integrals of the form 



After Feynman parametrization, change of variables b 2 — > m 2 — p 2 x(l — x), and 
k — > k — xp and Wick rotation we arrive at 






(2.2) 
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Figure 2.1: A next-to-tree- level Feynman graph of an arbitrary theory; it 
contains a loop and leads to divergent integrals. 



which is clearly divergent for large k (more specifically, this type of integral 
is called logarithmically divergent since a straightforward integration leads to 
logarithmic expressions) . 

At this point one can either condemn field theory altogether or think of a 
possibility to at least keep divergences out of observable quantities 1 . The first 
step of the latter must be to give a meaning to the divergent integrals. We will 
only be concerned with high momentum (ultraviolet) divergences since low energy 
divergences that appear if the loop contains a massless particle can be absorbed 
into soft-bremsstrahlung processes. Defining the integrals can be done in several 
ways, two of which are most commonly used today: Pauli Villars and dimensional 
regularization. For Pauli Villars regularization a cutoff A is introduced up to 
which the integral should be taken. This type of regularization is most often 
used in QED and we will therefore not go into any more details here. It shall 
suffice to add that the limit A — > oo has to yield a finite result for any observable. 
In QCD nowadays mostly dimensional regularization [21] is used. This works as 
follows: the integration of equation (|2.1j) is taken not in 4 dimensions but in D, 
where D = 4 — 2e. The value of the integral is then 



so that the divergence of the diagram is still present (as it has to be) in the 
pole of the T-function. Similar results hold for any other diagram, the pole in 
particular will always appear in the same way. In addition a mass scale \l is 
introduced, which is necessary to give the whole expression an integer dimension 
(alternatively, this also makes the coupling dimensionless) . fi can be interpreted 
as the scale at which the theory is defined. 



It is worth noting that singularities can even be encountered in classical electrodynamics 
of point particles. 
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(2.3) 
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2.1.2 Renormalization 



The basic idea of renormalization is astonishingly simple: just redefine parameters 
like masses or coupling as renormalized in such a way that the calculations yield 
finite results when carried out using the physical quantities. This is most easy to 
see on the level of the Lagrangian. Let the "original" Lagrangian of the theory 
be 

In terms of the "physical" parameters this Lagrangian will take the form: 

where "co" stands for the counterterms in the Lagrangian. It is tempting but 
untrue to say that one has "added" these counterterms. In fact, the Lagrangian 
has merely been separated into a divergent part and a renormalized finite one. 
The counterterms also correspond to Feynman diagrams, and their divergences 
will precisely cancel those of the original terms (the original terms now involve 
renormalized quantities, but the divergences remain, since nothing about the 
calculations changes). A theory is called renormalizable if a finite number of 
counterterms only is needed to render the results finite in all orders of the per- 
turbation series. 

Let us look once again at our test integral (|2.1|) . We have found that the 
result still includes the divergences as seen in the poles of T. We can expand the 
gamma function for small e and will, up to constant factors, arrive at 

(Any J dx(^j 2£ T(e) = ^ + \n4ir - j E - In ^ + 2 + 0(e) , (2.6) 
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in the massless limit, where 7^ is the Euler constant. This integral can for exam- 
ple be part of a two point correlation function, denoted by U(p 2 ). The correlation 
function can be split into a meaningless divergent part and a finite part that gov- 
erns the p 2 dependence. At this point it is necessary to impose a renormalization 
condition. While in QED on-shell renormalization is used in order to guaran- 
tee that the measured mass is the one that appears in the Lagrangian, this is 
less sensible in QCD, since quarks are not free particles and their masses cannot 
be measured. As for regularization there are multiple conventions, but the ones 
most often used are minimal subtraction (MS) and a modified minimal subtrac- 
tion scheme known as MS. In these schemes the correlation function is split in 
the following way, where we have again dropped irrelevant constants: 



r,-> 



n ren (p 



-JL 

ln47r-7ij-ln^ + 2 : MS 



- In =^r + 2 : MS 



_/ i + ln47r- 7£ ; : MS 

div ~ 1 1 : MS 



e 
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A third convention, that we will only mention, is the fi or Weinberg scheme 
where, by definition, the renormalized term contains only the logarithm. Note 
that the logarithmic dependence of the renormalized term is the same in the three 
schemes. All of the following work has been done in the MS scheme. 



2.1.3 The Renormalization Group 

The renormalization group equations can be seen as a connection between sta- 
tistical mechanics and quantum field theory. In short they govern the scaling 
behavior of the theory. This is dubbed "group" because of a groupoid structure 
|22[ ES] , even though it is not actually a group in the mathematical sense, since 
renormalization is not invertible. 

Now the basic statement is that all of the hard work we have put into the 
renormalization of the obtained expressions is highly artificial and should not 
be seen in measurable quantities. An observable should then neither depend on 
renormalization scheme nor scale. This statement of renormalization invariance 
leads to expressions that govern the scaling behaviour of physical quantities such 
as the correlation functions to be discussed in the next section. This, in turn, can 
sometimes also be true for non-physical quantities such as the Wilson coefficients 
of the operator product expansion, that will be the main topic of the next section. 
Consider, for example, a Greens' function denoted by T that shall depend on 
external momenta, the coupling constant and masses. The relation between the 
the bare Greens function and the renormalized one is then 

T({pi}, a, m a , n) = lim Z^fj,, e)T Q ({pi}, a , m afi , e) , (2.7) 

where Z^/j,, e) is the appropriate combination of renormalization constants. It is 
customary to define a = a s /Tc and we will do so throughout this work 2 . Taking 
the total derivative of Eq. ()2.7|) with respect to /i and keeping in mind that the 
bare Greens function must definitely not depend on /i we obtain 



Q 

P a zZ ^lo,m a - + 7 Z 



T({pi},a,m a ,n) = 0. (2.1 



This is known as the Callan-Symanzyk or renormalization group equation and 
we have introduced the beta function and anomalous dimensions as 

da fi dm a /i dZ 

P = 3" ' la = J- , lz = ~ -i- ■ 2.9 

a/i m a djj, Z djj, 



Let us first discuss the running of the coupling constant. In the MS scheme the 
(3 function depends on a only (it is therefore called a mass independent scheme) 

2 In general a Greens function will also depend on the gauge parameter, unfortunately often 
also denoted by a; we will use Landau gauge and thereby omit this factor. 
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and can be expanded 



/3(a) = fca' + foa 6 + fca? + ... , (2.10) 

where nowadays all the coefficients up to (3± are known. (3\ and fa are, in fact, 
scheme independent. All known coefficients are given in the appendix. An explicit 
expression for a can be obtained by integrating the defining equation. Up to 
leading order the resulting expression reads 

a (Q) = 14 .fl a M, a - (2 - n) 

1 + /3 1 o(/x)lnif 

In order to remove the renormalization scale one conventionally introduces a 
another mass scale A in such a way that the expression simplifies to 

a(Q2)= fti^r (2 ' 12) 

For higher orders a further expansion in powers of logarithms holds. The scale 
A is at several hundred MeV and will depend on the number of quarks included 
in the theory which causes it to change at quark thresholds. This is achieved 
by imposing matching conditions for the effective theories with the appropriate 
number of quarks. Conventionally a is given at the scale of the Z boson mass for 
which the Particle Data Group value is 



a(M z ) = 0.119 ±0.002. (2.13) 

From Eq. 1)2.12)1 we can see that for positive (3\ (as is the case for "real life QCD" , 
see appendix A) the coupling exhibits both confinement and asymptotic freedom 
and the scale A becomes meaningful as the scale at which confinement is effective. 

In this work we will define the coupling implicitely, by directly integrating 
Eq. (|229). We arrive at 

/ v a(Q) 

nfHw< (2i4) 

a(fi) 

which can be numerically solved for given values of the /3-function coefficients. 
We will be using the known 4-loop expression that should be a very good approx- 
imation, unless the j3 function has some unexpected higher-order pathologies. 

The treatment of the running mass is in perfect analogy to that of the running 
coupling: again the anomalous dimension 7 can be expanded into a power series 
in a 

H dm 

7(a) = — = 7ia + 7 2 a + ... , (2.15) 

m dfi 
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where the coefficients are known up to fourth order and given in the appendix. 
A simple integration of (j2.15|) gives a series for m in powers of a. Up to next-to- 
leading order this is: 



m(Q) = m(/i) 



HQ)' 


71 /ft ^ 


a(/j,)_ 





1 + HQ) - a(n)] (j^ - ^rj^j • (2-16) 



The coefficient fix is positive as noted above and the same holds for jx, so we 
conclude that masses will decrease along with the coupling constant for high 
energies. 

As for the coupling we will choose to solve the resulting integral equation 
numerically. It is 

tn\ a(Q) i \ 

^ dx p (2 . 17) 

<haO 

Taking the exponential of both sides of this equation shows that the ratio of two 
quark masses at the same scale is scale- independent: 

a (Q) . . 

m{Q) = m(fi) exp f dx^\ , (2.18) 
J p[a) 



a(p) 



since the exponential factor on the right hand side will drop out if we take the 
mass ratio. 

For the light quarks (u, d, s) it has become customary to quote values at either 
2 or 1 GeV, since perturbation theory breaks down at about this scale. On the 
other hand, masses for heavy quarks are usually given at the scale of their own 
mass. The current values for the quark masses are then |2~H] : 



m u (2GeV) = 2.9±0.6MeV, m d (2GeV) = 5.2±0.9MeV, 
m s (2GeV) = 98±18MeV, m c (m c ) = 1.28 ± 0.09 GeV, 

m b (m b ) = 4.21 ±0.06 GeV, m t {m t ) = 165 ±5 GeV. 

2.2 Operator Product Expansion (OPE) 

The operator product expansion was introduced in 1969 by Wilson j2Sj as an 
alternative to the current algebra methods used at the time. He proposed to ex- 
pand nonlocal operator products into a sum of local operators, thereby relegating 
the nonlocality into the expansion coefficients and thus separating long and short 
distance effects: 

A(x)B(y) x ~ y £ C n (x - y)O n (y) , (2.19) 
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where A, B and O n are operators and the C n (x — y) are known as the Wilson 
coefficients of the OPE. For simplicity y is usually set to zero, a convention that 
we will adopt from here. It is, however, important to note that the above equality 
is valid only in a weak sense, meaning that the operators have to be sandwiched 
into initial and final states. In general the sum on the right hand side involves 
operators of different mass dimensions, which causes the coefficient dimensions to 
vary accordingly. As we will see it turns out that, in the specific cases discussed 
in sum rules, this leads to a suppression of higher dimensional operators, a fact 
without which the OPE would not be of much use. 

Since the product of two operators at arbitrarily small distances can in general 
be singular one expects the same to be true for the coefficients. However, the 
nature of the these singularities is governed by the underlying exact and broken 
symmetries, namely scale invariance; for a complete discussion the interested 
reader is referred to Ref . or [23] . At this point we will only mention that the 
scaling dimension (as opposed to mass dimension) of the coefficients is 

C n (x) x~ x ", A n = d{A) + d{B) - d(O n ) , (2.20) 

just as one would expect from dimensional analysis, only that d(X) is not the 
conventional mass dimension d c but 

d(X)=d c + 1 (a*), (2.21) 

where 7(a) is the anomalous mass dimension of the operator and a* is the fixed 
point of the renormalization group equation. Operators that generate symmetries 
as well as the identity operator 1 among others will keep their naive dimensions 
through this process. 

The operator product investigated in sum rules is the two point correlation 
function of two-quark currents: 

U^(q) = 1 J d 4 xe^(Q\T(f(x)f(0))\Q) . (2.22) 

Here the current is built from any two quark- ant iquark flavors. 

f(x) = q i Tq j (x), (2.23) 

with T any combination of Dirac matrices introduced in order to obtain the 
desired quantum numbers. The coefficients for various currents have been calcu- 
lated and can be found in the literature (see e.g. Ref. ^0] for a compilation of 
coefficients and references to the original works). Taking for example 

f(x) = - (: u^u{x) - d^d{x) :) , (2.24) 

allows one to model the p-meson with rather good accuracy. As usual, the no- 
tation : : implies a normal ordering. Indeed, this vector current sum rule is 
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the standard example of an OPE and has already been discussed in the original 
paper by SVZ [H1IZ|, while the explicit calculation to this correlation function is 
presented in great detail in Ref. [211] • The lowest dimensional scalar operators 
(Dim < 6) that will appear in the OPE are 

1, mm, F^F^, m t qa^F^ a q, qVqqTq. (2.25) 

A normal ordering is understood in all of these operators. In ordinary perturba- 
tion theory these vacuum expectation values vanish by definition except for the 
unit operator. For the others they are introduced as "condensates" in nonper- 
turbative QCD in order to describe the nontrivial vacuum structure. A natural 
origin of the quark condensate can also be found in chiral perturbation theory, 
where it arises as the order parameter of spontaneous chiral symmetry breaking 
and is then responsible for the pion mass (the famous Gell-Mann Oakes Renner 
(GMOR) relation) [2Z1I2H1- At this point an additional problem turns up: the 
standard proof of OPE relies on Feynman diagrams and perturbation theory, so 
its use in a nonperturbative regime is far from obvious. This causes the expansion 
series to break down at a certain dimension d cr a that has already been estimated 
by SVZ. In practise the calculations are usually performed up to dimensions 6-8, 
depending on the strength of the suppression. 

Now let us see how the higher dimensional operators arise in a semiperturba- 
tive expansion. First, however, we need to decompose the correlator into parts 
of different Lorentz structure by way of the Ward identities. This can be done as 
follows: 

W\q) = (gY-^n^ + ^^T^^V) (2-26) 
„„ (mi t mi) r / v , m 

+ q2 [{QiQi) T (gjgj)] ) 

IT ' p (q 2 ) being the correlator of scalar (pseudoscalar) currents. The proof for this 
relation can be found in Ref. ;29j but the structure can also be understood intu- 
itively: the term in the first brackets is the "ordinary" vector current correlator 
that implies transversality, the second can be understood by looking at the corre- 
lator of the divergence of a vector current. Using the Dirac equation this can be 
transformed into a scalar correlator multiplied precisely by the given expression 
in the masses. The last two terms are subtraction constants for the scalar current. 
All of this will be discussed in more detail when we study the scalar two point 
function in chapter 3. In our example of the vector current from above only the 
first term contributes, and we will begin by sketching the calculation for the unit 
operator (Fig l2.2l a)). Contracting both sides of Eq. (j2.26j) with g^ u we get: 

U(q 2 ) = — ? / d D xe iqx (Q\T(: u(x)^u(x) - d(x)^d(x) : (2.27) 

Aq 2 (d — 1) J 

x : m(0) 7m u(0) - rf(0) 7M rf(0) 
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q-p 




Figure 2.2: Quark condensates: pictorially the nonlocal condensate b) corre- 
sponds to cutting one propagator in the Feynman graph for the perturbative 
operators a). 



to which we apply Wick's theorem: 

1ST j f fj^T 

n (? 2 ) = - 4g 2 (D C _ 1) J J^nlTr[YS u (p) lfl S u ( q -p)]+u^d}. (2.28) 

Completing the calculation gives a coefficient of gf^d — ln(— j^)) where /x is the 
renormalization scale that arises in dimensional regularization. 

But let us go back one step to Eq. ()2.27|) and leave one quark- ant iquark 
pair uncontracted. Pictorially this corresponds to cutting one propagator in the 
Feynman graph (see Fig. 12.2b )). One is left with expressions including terms like 
(:q(x)q(0) :) that can be Taylor expanded in order to arrive at well defined local 
condensates (:g(0)g(0) :), often written as (qq). Note that this simple derivation 
naturally leads to normal ordered condensates, a subtility that becomes important 
in calculations: writing the OPE in terms of operators that are renormalized 
in a non minimal scheme (i.e. normal ordering) leads to coefficients including 
powers of mln(^-). The usual procedure is to set — q 2 = /i 2 because all of the 

2 

logarithms of the type ln(^-) are then set to zero (this procedure is known 
as "summing up" the logarithms and can be done because, as we will see, the 
correlator is related to an observable quantity and has to obey a renormalization 
group equation). Unfortunately, this summing up process does not get rid of 
mass logarithms and in fact causes them to be large, which makes calculations 
of nonleading contributions impossible. Mass logarithms are then interpreted 
as remnants of long distance effects so that the separation of long and short 
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distance effects in non-minimal schemes is incomplete and working in minimal 
schemes appears advantageous. A more elaborate discussion of this subject can 
be found in Ref. [201; where the relation between normal ordered and minimally 
subtracted condensates are also given. 



2.3 Dispersion Relations and Sum Rules 

This section will be dealing with the analytic properties of the two point cor- 
relators introduced above and their consequences. It is shown for example in 
Ref. [31 J that a general two point function can be written as an integral over a 
spectral function p(s) 

oo 

U(s) = [ dtp(t) _. , (2.29) 



where the spectral function is the sum over all possible final states A: 

p(q) = (2vr) 3 ^ \(n\j x (0)\\)\ 2 8*(q-p x ). (2.30) 

A 

We have also introduced the notation q 2 = s that we will be using frequently from 
now on. The correlator is then an analytic function except for possible poles on 
the positive real axis. Now remember that according to Cauchy's theorem every 
analytic function f(z), that we will also assume to fall off fast enough for s — > oo, 
obeys the following relation for any closed curve C with t as an interior point: 

f(t) = — [ dz^- . (2.31) 
Jy ' 2iriJ z-t y ' 

c 

Let C be the curve from Fig. 12.31 it consists of a circle of radius R around the 
origin, with an interval of 2s left out around the real axis. The curve runs above 
and below the axis and is finally closed by a circle of radius s around the origin. 
Integrating the integrand of Eq. ()2.29|) over this curve and applying Cauchy's 
theorem, we find: 

dt — P ^ . = 2mp(s + is) (2.32) 



■is 



dt + / dt ^ + f dt f (f) 

t — s — is J t — s — is J t — s — is 

?i C2 C3 

+ / *■ p(t) 



t — s — is 



The integral over the large circle vanishes for R— > 00, as does the integral over 
the small one, if the integrand is analytic into the origin, which we will assume. 
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Figure 2.3: Integration contour for dispersion relations. 



The remaining terms are 



2mp(s) = 2 J dtp{t)m8{t- s) (2.33) 
o 



lim / dt 



e-*0 J t — S — IS' 

C1+C2 



or, written in a more symbolic form: 



' V \- m8{t - s). (2.34) 



t — s — is t — s 

where V indicates that the principle value is to be taken. We can now apply 
Eq. ()2.34j) to the correlator and compare real and imaginary parts of both sides. 
Beginning with the imaginary part, we find 

p(s) = -Imn(s) . (2.35) 

7T 
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Analogously, the real part is given by 



ReII(s) = f dt-^-, (2.36) 
o 

where convergence of the integral will have to be checked. The same holds for 
the integral in equation Eq. ()2.29j) . which may not converge either if p(s) does 
not fall of fast enough. In order to make the integrals convergent one can expand 
both sides into a Taylor series and subtract the first term before integrating (this 
is generally referred to as "subtractions" ) : 

oo 

11(a) - n(so) = / dtp® (j^— s - j-j-} (2.37) 



oo 

(s — s ) J dt 





s , 

Pit) 



(t-s)(t-s ) 



This integral will obviously converge better due to the higher power of t in the 
numerator. In certain situations more than one subtraction is needed (the number 
of necessary subtractions can easily seen by power counting since the integrand 
needs a negative mass unit to converge). The generalization of the above result 
can then be shown by induction to be 

U(s) = g a k (s - s ) k + (s- so) n J dt , (2.38) 



fc=0 



where 



1 d k U(s) 



, , , , ■ (2-39) 

k\ ds K s=s o 

A very simple sum rule can already be obtained by taking the OPE of the cor- 
relator on the left side of the equation and using the measured spectral function 
on the right side. Often the spectral function can be directly related to a cross 
section or decay rate by the optical theorem: in the case of the ratio of quark 
anti-quark pair production to muon pair production in e + e~ collision for example 
this relation is given by 

#( s ) : = ae+e ~^" = 127rImII(s) . (2.40) 

® e+ e~ — > /U+ fi~ 

At this point it is worthwhile to pause and recapitulate what we have essentially 
done: we have, so to speak, built a bridge from phenomenological parameters (the 
spectral function) to the fundamental parameters of QCD (coupling constants 
and quark masses). Vacuum effects have been included in a rather small number 
of parameters that represent static vacuum properties. It is crucial to assume 
universality for these parameters in order to make sensible predictions. 
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Applying this sum rule to our "toy model" of the p meson turns out to be 
a bit disappointing because of conflicting choices for s on each side. We expect 
the spectral function to be dominated by the p meson at low energies and would 
thus like to choose a very small s, while on the other hand the OPE involves 
terms that are suppressed by factors of quark masses over momenta which would 
favor the choice of large s. The first idea would be to use a compromise value 
of, say 1 GeV 2 , which is okay for the OPE but still not good enough as p meson 
saturation approximation. Clearly the sum rule will have to be improved. There 
are different ways of doing this and we will discuss them in the next section. 



2.4 Different Types of Sum Rules 

We have seen above that the simplest sum rules can often be only mediocre 
approximations and we will now try to improve on them. There are basically 
three types of sum rules that are widely used |HT]- They all differ in the way 
the analytic properties of the spectral function are exploited or certain limits are 
taken. 



2.4.1 Borel Sum Rules 

This is the kind of sum rule that was used for light quarks in the original paper 
by SVZ and enhances the low energy region of the spectrum. We will, however, 
choose another motivation: going back to the vector meson sum rule of the last 
section for a minute, dimensional counting shows us that we need one subtrac- 
tion constant to arrive at a well defined expression. In a more general case the 
correlation function can have dimension d and the asymptotic behavior is 

lim -Imn(s) => As d (1 + ca s {s)...) . (2.41) 

s — *-oo 

Then d+1 subtraction constants are needed. The subtraction constants generally 
depend on the renormalization scale and drop out of any physical observable, but 
for calculations it is necessary to get rid of them. The easiest way to achieve 
this is to take d + 1 derivatives with respect to Q 2 = —s. We are, of course, not 
limited to d+1 derivatives, since any number N > d+1 of derivatives removes the 
subtraction constants. The resulting expression is sometimes called the moment 
sum rule: 



{N-d-l)\^ ' (dQ 2 ) N ' 



o 



(2.42) 

{t + Q 2 ) d+1 \t + Q 2 ) 7i 



, N(N-l)---(N-d) ( Q 2 \ {N - d) i T . . 
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To get rid of any arbitrary number of subtraction constants we can take the limit 
iV — ► oo; additionally we will let the energy Q 2 go to infinity so that the ratio 
= M remains fixed. This corresponds to taking the Borel transformation, 
defined as 

1 / d \ N 

B -J}^J^fW) { - Qr W) • (243) 

of the first well defined moment. It can be shown (see Ref. and references 
therein) that the Borel transformation is the inversion of the Laplace transform 3 
so that the Borel transformed dispersion relation reads 

oo 

M 2 B U(Q 2 ) = J dt p{t) e~ t/M2 . (2.44) 
o 

It is easy to see that the low energy region of the spectral density is enhanced, 
as mentioned above, so that this type of sum rule is particularly suitable when 
studying low energy resonances. In practise the phenomenological spectral func- 
tion is known only up to some threshold energy sq. Above this threshold the 
imaginary part of the OPE for the correlator is taken as a theoretical spectral 
function. If the threshold energy is large enough this should be a rather good 
approximation. There have already been several successful attempts to extract 
light quark masses from Borel sum rules; the results for the strange mass shall 
be reviewed and discussed at the end of this work. 



2.4.2 Finite Energy Sum Rules 

For finite energy sum rules (FESR) consider again the contour from Fig. 12 .3[ only 
this time we will not let the radius of the large circle go to infinity. Instead, we 
will integrate over a circle of finite energy so (hence the name). The function to 
be integrated is a weighted correlation function: 

dsw{s)U{s), (2.45) 



2vri 
c 



where w(s) is an arbitrary analytic function. According to Cauchy's theorem the 
integral over the whole contour has to vanish. The parts of the contour along 
the axis will pick up the discontinuity in the spectral function (p has to obey 
Schwarz' reflection principle) and we find, if we split the contour into two parts: 

I dsw{s)U{s) = J dsw{s)p{s). (2.46) 



2ni 

s=sq 



3 This fact leads to the name Laplace transform sum rules that can also be found in the 
literature instead of Borel sum rules. 
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The obvious advantage of this approach is that the spectral functions needs to be 
known only up to the energy y/so, while for Borel sum rules it has to be known to 
arbitrary high energies. Unfortunately, the convergence behaviour of FESR may 
not be as good as that of BSR, depending on the weight function. The question 
of subtractions in finite energy sum rules is dealt with by partial integration: 

j> dsw(s)U(s) = - j> ds[W(s)-W(s )]^^-, (2.47) 

s=so s=so 

where W(s) is the integral over the weight function. For higher order subtraction 
polynomials this procedure can be repeated and all subtraction constants will 
then be successively removed. Further subtleties of FESR will be discussed in 
the next chapter that shall deal exclusively with the strange mass FESR. 



2.4.3 Gaussian Transform Sum Rules 

Here we will consider jH2] the Gaussian transform of the spectral function 

G(s ' T)= v&/ dfexp r^r (i) - (2 - 48) 

Depending on the choice of s and r any energy region can be probed, so one 
is free to investigate the region that is most easy to parametrize. In order to 
calculate G we will define the following combination, where we assume that the 
dispersion relation has at most one subtraction constant (the calculation can 
easily be generalized to arbitrary subtraction polynomials): 

"(*) = m (2 49) 

(t-s) 2 + A 2 (t-s + iA)(t-s-iA)' y ' 1 

Once again this function is integrated over our contour from Fig |2.3l with R — > oo, 
and we find 

nt-MA) , u { s-iA) _ J >n W 



lA -iA J (t-s) 2 + A 2 

v ' 

Now we can apply the Borel transformation with respect to A to obtain the final 
result: 

i n A II(s + zA) n(s-iA) 1 7, / (s-t) 2 \ 
,2tB + y ' = -== / dt exp -i '— p(t) 



iA -%A y/l 7tt J V 4r 

= G{s,t). (2.51) 



Looking once again at the definition of the Gauss transform we observe two things: 
the first is that for r equal to zero the Gaussian transform merely reproduces the 
spectral function, the second is that it obeys a differential equation: 

d 2 G(s,r) _ dG(s,T) 
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This equation is known in another context as heat evolution equation and we can 
try to apply this analogy to our spectral function. Taking r as our analog to a 
time variable, we can call the spectral function the initial condition. In pQCD the 
Gaussian transform is an expansion in a(y/r) which is not calculable for small r 
because of confinement and thus comparisons of hadronic spectral functions with 
pQCD cannot be made. But one can evolve the spectral function "in time" by 
going to larger r and then compare the evolved spectral function with predictions 
from the corresponding pQCD expressions. 
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Chapter 3 
r-Decay Sum Rules 



Hopefully you will understand the simple examples, be- 
cause if you do you will understand the generalities at 
once - that's the way I understand things anyway. 
(Richard P. Feynman, Elementary Particles and the 
Laws of Nature) 



3.1 r-Decays and Finite Energy Sum Rules 

Finite energy sum rules were first applied to r-decays systematically by Braaten 
et. al. |HSj when studying inclusive r-decays, following the ideas of Schilcher and 
Tran [HI]. Since their work is the basis of our studies and we need to fix our 
notation, we will review it briefly, sometimes taking a more modern point of 
view. At the same time the problems of this approach will become evident. It is 
interesting to note that at the time these sum rules were used to determine the 
QCD coupling at the scale of M 2 , another example of how the application of sum 
rules often goes beyond the original intention. 

Throughout this section we will be considering the ratio of hadronic to leptonic 
r-decays defined as 

Rr = Ti Z ~ hadr ° m) E ^ 3-647 ± 0.014[IH! . (3.1) 
A first estimate gives, neglecting strong and elektroweak corrections: 



Rr = N, 



V ud \ 2 +\V US \ 2 ~3. (3.2) 



Because of the inclusive nature of r-decays this ratio can be decomposed into 
several contributions associated with different quark currents. These are vector, 
axial vector and strange contributions. The strange decay ratio cannot be de- 
composed further because parity is not a conserved quantum number in weak 
interactions. We can now write the Lorentz decomposition (|2.2fij) in a slightly 
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different way, where we will also absorb the subtraction constant into the scalar 
term: 

= {q»q v - q 2 g^)rt T \q 2 ) + gVTI (L) (g 2 ) . (3.3) 

The relation between the component Ii^ L \q 2 ) and the scalar spectral function 
IlW)(g2) from Eq. §F2B 

can be obtained by multiplying both equations with 
q^q u . For Eq. ()2.26|) we find: 

q^q v UT(q) = K T m j ) 2 n 5 ' p (g 2 ) + (m, T m 3 )[{q iqi ) T (m)] , (3.4) 

and for Eq. Q : 

g^n^(g) = n( i )(g 2 ) S 2 , (3.5) 
so that finally the relation reads 

s 2 rt L \q 2 ) = { mi =F m,) 2 n 5 ' p (g 2 ) + (m, =f m^ffe) =F (fc^)] • (3.6) 

The second term on the right hand side is, as mentioned before, merely a sub- 
traction constant, but the first term will be needed later for the relation between 
the spectral functions. Additionally, we need the relation between the transversal 
and the vector component. This can be seen by adding and subtracting q 2 g fJ,u in 
equation ()3.3|) : the result of this exercise is 



U v ' A (s) 



(3.7) 



Using the optical theorem 1 the decay rate can be written as an integral of the 
imaginary part of this correlator 



J/r 



R T = 12tt 



ds I s \ 
M 2 \ ~M?) 



| ^i + 2^imn( T )( s ) + imn( i )( 



(3i 



This expression looks very much like the right hand side of Eq. ()2.46|) . but before 
we go on we have to discuss weight functions in FESR. The problem is that the 
integral runs along the real axis where perturbation theory completely breaks 
down and the OPE is no longer valid. Since we will want to turn the above 
expression into one that involves an integral of Il(s) over a circle around the origin 
we had better make sure that the contributions along the axis are suppressed. 
In case of the r-decay we naturally arrive at a weight function that meets this 
requirement, so we could have proceeded without thinking about this topic at all. 
It is, however, necessary to know what one is doing when working with general 
FESR. That said, we can now go on along the reasoning that led to Eq. ()2.46|) 
and obtain: 

2 r / \ 



R T = 6ni 




|s|=M| 



M 2 ) K 1+ V 2 , 



n( T )( s ) + n^( s ) 



(3.9) 



lr The reader is invited to convince himself of the correctness of this formula by following the 
proof given in appendix |0 
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We see that the integral goes over the full circle but suppresses contributions 
along the axis by a double zero. Experience shows that this is enough to make 
the sum rule well satisfied. Sum rules of this type are called pinch weighted finite 
energy sum rules or pFESR. 

The vector snd scalar contributions can be cleanly separated by rearranging 
the terms: 



R T = Qni j 



|s|=M| 



ds ( s 



(3.10) 

To get rid of renormalization scheme dependent subtraction constants 2 we define 

B< L+r >l» = -s^(n< L+T >(s)) , (3.11) 

and use partial integration to arrive at 



|s|=Af? 



(3.12) 



which can be calculated in QCD. 

Further information on different parts of the spectrum can be obtained by 
defining moments [SHI, that include additional weighting factors: 



s \ I s \ dR T 



R " = I \Mj) ~da dS ' (313) 

Just as above these moments can be converted into weighted contour integrals 
over the derivatives defined in ([3.11)1 : 



Ru =-*i f ^ (3^ T (^ (i+T) (^)+4^W^(^)) , (3.14) 
1*1=1 



where the functions JF fe '(a;) contain the weight functions as well as all the factors 
that arise because of partial integration. It turns out that for nonzero values of 
I the sum rule is not very well satisfied, so we will study only I = moments for 
different k. The relevant functions J rkl (x) are given in the appendix as well as a 
general representation. 

At this point we need to remember that we have decomposed our decay rate 
into several contributions, as well as to keep in mind that the strange channel 

2 Dimensional analysis shows that one subtraction constant is indeed enough. 
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is Cabibbo-suppressed, while the axial and vector channels are allowed. We will 
also include the conventional electroweak radiative correction terms for which the 
numerical values is Sew — 1-0201 ± 0.0003 [36J |37J [SH]- The correction term 5' EW 
[3*7] is absorbed into Sew- 

Inserting the operator product expansion we then obtain 



Rr = 3 [\V ud \ 2 + \V US \ 2 ] Sew\1 + 5 (0) + E H 2 d cC ] + sm 2 c ti£>] 

D>2 

(3-15) 

where 8 C is the Cabibbo angle. The coefficients of the operator product expansion 
are "hidden" in the summed terms d\f , which also include the integrals over 
weight functions and so on. We will be concerned with the precise structure 
of these terms in the next section of this chapter. Even now we can notice 
several details: the dimension-zero term does not involve quark masses and is 
then independent of the flavors. It only "sees" them through the elements of the 
CKM matrix (or in fact is does not because the sum of the two squared elements 
is very close to one). The perturbative contribution as well as all other SU(3) 
invariant effects such as instantons will drop out, if we only consider the SU(3) 
breaking part of R T : 

SRr = - ^ = SS EW £ [6™ _ 6S)] > (316) 

I 'ud I I *us | D>2 



which then allows us to considerably reduce the errors. 



3.2 The OPE for SU(3)-breaking Contributions 

In this section we will collect all the terms needed for the OPE side of the sum 
rule. We will only give general expressions and refer to the appendix for the 
explicit coefficients. Additionally, the dimension-zero perturbative contributions 
are omitted, since they cancel if we study SU(3) breaking effects. The interested 
reader is referred to Ref. |3"I] for a detailed discussion of these perturbative con- 
tributions. The only dimension-two operators that can be constructed are masses 
so that the dimension-two corrections to the power series will give the needed sen- 
sitivity to the strange mass. We will separately present the longitudinal and the 
transversal components of the spectral function, or equivalently, as we have seen, 
the scalar and the vector components. Especially, we will need to discuss the 
scalar contributions very carefully, for reasons that will become apparent later. 
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3.2.1 Dimension-two Corrections for the Scalar Contribu- 
tion 



The dimension-two terms could in principle consist of any combination of squared 
mass terms so it is instructive to first consider the form that the resulting terms 
are expected to have. To do that we need to discuss the correlator of scalar 
currents; it is more convenient, however, to start from the correlator of diver- 
gences of vector and axial vector currents respectively, since we are then working 
with renormalization invariant objects. To see how they are related to scalar 
correlators we will show more explicitely what we have already mentioned below 
Eq. ()2.26|) : applying the Dirac equation to 

j'(x) = d»(s lfl u{x)) , f(x) = d»(sw 5 u(x) (3.17) 

immediately gives 

f(x) = i(m s - m u )(su(x)) , f(x) = i(m s + m u )(s^ b u(x)) , (3.18) 

so that the correlator should appear multiplied by the sum or the difference of the 
masses squared (we could have gained this insight by simply looking at Eq. (j3.6|) 
but now we understand better how the scalar correlator appears as part of a 
vector one). The dimension- two contribution is then 



D^\s) = -^ W2 (m l Tm ] ) 2 Y,dn(q 2 /^)^\ (3-19) 



where a and m are taken at the renormalization scale fi. The coefficients d n (q 2 / /x 2 ) 
are in general a sum over powers of logarithms: 

n 

d n (q 2 /lx 2 ) = Y,dnjV , (3.20) 

3=0 

where L = ln(— q 2 /jj?). The quantity D^ L \s) is related to an observable (the 
spectral function), so it has to obey a homogeneous renormalization group equa- 
tion, which leads to relations between the d n j. The RGE is: 

(-4 - - -< {a) ( m -ik + ) °? w = ° • < 3 ' 2i > 

where the factor of two comes from the fact that L involves fi 2 . Inserting 
Eq. ()3.19j) with the representation ()3.2())1 into the RGE and comparing the terms 
of same order in L and a gives the following scaling relations: 

du = -jid 00 , (3.22) 

1 r 



d 



21 



(#l + 271) doo + 27 2 d. 



00 
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d 



22 



-71 (A + 271) d 00 



tZ: 



1 r 



31 



2 (ft + 71 ) d 20 + (& + 2 72 ) d 10 + 273d 



00 



32 



d 



33 



(I1P2 + 272(271 + Pi)) d 00 + (Pi + 71) (27i + Pi)d w 
~li (A + 7i)(27i+/3i)d 00 - 



This means that all that needs to be calculated are the d,o coefficients; all others 
are completely determined by the scaling laws of the renormalization group. We 

->L,2 



could now go on and just integrate the above expression for D^'y/^s) over our 
circular contour. Then we have to set the renormalization scale at fi 2 = M 2 . It 
turns out, however, that this leads to a rather bad convergence behavior due to 
large logarithms when integrating over the circle. This can be avoided by using 
the so called "contour improved perturbation theory" (CIPT), which means the 
following: both the renormalization group and the operation of integrating are 
linear, and we can set our renormalization scale as fi 2 = —q 2 before doing the 
integration. We will thus sum up the contributions point by point along the 
contour. As the masses and couplings are given at the renormalization scale the 
resulting expression will now involve integrals over both the running coupling 
as well as the running mass. This not only leads to simpler expressions for the 
integrand but also to a better convergence, as has been shown in Refs. jHUEO]. 
For illustration purposes we will quote here the perturbation series as given in 
Ref. 110] for ordinary perturbation theory and CIPT. The quoted values are the 
sum of the longitudinal and the vector part that we have not yet discussed. The 
numerical values are: 



1 + 0.5866 + 0.5566 + 0.5150 + ... 
for ordinary perturbation theory and 

0.9840 + 0.4613 + 0.3421 + 0.298 + ... 



(3.23) 
(3.24) 



for CIPT. The better convergence of the latter should be obvious. 

Finally we can write 5R T in a form that is more transparent for the numerical 
evaluation: 

m 2 (M 2 ) 



SR 



L.2 

kl 



6S 



EW 



n=0 



L(n) 
kl 



where 



B 



L(n) 
kl 



2m J x Ly ' 

|x|=l 



m(-i 2 M 2 T x 



m{M 2 



a n (-^M 2 x 



(3.25) 



(3.26) 



and Ed is the ratio of of strange and down quark masses calculated rather accu- 
rately in Chiral Perturbation Theory [S]: £d = rrid/rns = 0.053 ± 0.002. Later we 
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will also encounter e u , for which the value is e u = m u /m s = 0.029 ± 0.003. We 
have also introduced a parameter £ that is of order one and will be varied later 
as an estimate for errors. This is necessary because the truncation of the per- 
turbative OPE contributions leaves a residual scale dependence in the obtained 
expressions (the reason being that the RGE for the correlation function "mixes" 
different orders of a and L, as can be seen from the scaling relation for the Wil- 
son coefficients). In the language of s and /i we find that £ 2 = — fi 2 /s so that all 
logarithms are zero for £ = 1. 



3.2.2 Dimension- four Corrections and Higher 

While there was only one operator of dimension two, namely the mass, we will 
have to consider several contributions with dimension four. First of all there 
is obviously a term with the fourth power of the masses, next we have to take 
into account quark condensates multiplied by a mass and finally there could 
also be a gluon condensate operator. Again we can begin by trying to get an 
intuitive feeling about the structure that we ought to expect. Remember that 
the longitudinal contribution is just the scalar correlator times a squared term 
in the masses so that the dimension-4 terms of the longitudinal contribution will 
be related to the dimension-2 contributions of the scalar series. This means that 
there are no condensate terms, except for those contributions from the subtraction 
constants in Eq. (|3.fij) . A general expression for the dimension-2 series of the scalar 
correlator is given in Ref. p|, from which we can easily calculate the following 
expressions relevant to us: 



D LA 



M 2 s 
3 



((mi T nij) (qiqi T qjQj)) (3.27) 



82 112 



8tt 2 

x (4 + 4 m£ + a(-£ 2 s) ( - 8C(3) + y + ^ In £ + 16 In 2 £ 
±m i (-e 2 S )m,(-£ 2 S )(6 + 41n£ + a(-£ 2 s)(48 - 8((3) 
+^ln£ + 161n 2 

that has to be integrated over the contour according to Eq. ()3.12j) . While the 
masses will be scaled according to the renormalization group equations for this 
integration, we will take the combination of condensate times mass to be scale 
independent. This is because they are given to a very good approximation in 
chiral perturbation theory by the GMOR relation and should then indeed be 
scale invariant. 

Next we have to include terms of dimension 6. The same rationale that fixed 
the structure of the dimension 4 can again be applied and shows that we have 
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to take into account every possible operator of dimension 4 and multiply them 
by the standard mass prefactor. We will then consider m 6 terms, condensates 
multiplied by cubic mass terms and finally a gluon condensate term. The relevant 
expressions, calculated again from those given in Ref. are: 

3 (m^-es) T mft-fsj) ( ± m, 3 (-£ 2 S K(-£ 2 S )(3.28) 
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Finally we will include operators of dimension 8. The complete expressions will 
obviously become quite lengthy, but we will only consider the terms that are 
numerically most important. The largest would be a leading order 4 quark con- 
densate but this turns out to be zero jH|, so we shall include the mixed condensate, 
that is suppressed by only one mass power (apart from, obviously, the squared 
terms in the prefactors) and the first order four quark operators. Since these four 
quark operators cannot be reliably estimated, they are conventionally reduced to 
the chiral condensates using the vacuum saturation approximation 

For that we shall take a general example of a four quark condensate with 
arbitrary Dirac structure T and fixed color structure t a , where t a are the Gell- 
Mann matrices that generate the SU(3) multiplied by an overall factor of one 
half, in other words we discuss 



(3.31) 



or 



^l,iy^ 2 p\^kl^mn(^k^ul^pm ( lXn) • (3.32) 

Dirac indices are denoted by Greek letters and run from one to four, while the 
Roman indices, corresponding to the color, run only from one to three. A sum 
over flavors, denoted by capital Roman indices, is not included. 

We will now insert a complete set of eigenstates into the bracket and assume 
that the contribution from the vacuum state is the largest. Since the vacuum is 
a flavor- and colorless scalar state, we arrive at 



■pl T-.2 .a .a x °HV°kl x & pX &mn , A „B\ / -C 
1 pu l pX t kl t mn°AB— ] ~®CD ' " 



4 3 



4 3 
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(3.33) 
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where the factors in the denominator have to be included to ensure proper vacuum 
normalization. This expression can be more compactly written as: 

'" Trr 1 Trr 2 Trt a Trt a 5 AB (5cD(g A g B )(g C g D ). (3.34) 
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The trace of the color matrices is zero, but we will keep track of the explicit ex- 
pressions, since the final result we will obtain can then more easily be generalized. 
Additionally, we have to consider the expression where the quarks in ()3.33j) have 
been exchanged, i.e. 

— I" ' p \tkl^mn('l^k ( lXn ( lpm ( lul) 5 (3.35) 

where the additional minus sign stems from the interchange of the fermion fields. 
Again we insert the vacuum and find 



- r ^r 2 A t^C n 5AD^^5cB^^(g^gf n )(^^) > (3.36) 



or 

- ^tV x V 2 TYt a t a 5 AD 5cB(q A q D )(f^) , (3.37) 
Combining the two we eventually arrive at 

^{TvY 1 TiY 2r Tit ar Iit a 5 AB 5cD{q A q B ){fq D ) (3.38) 
-TrV l V 2 lktH a 8 AD 8 CB {q A q D ){fq D ))- 



Finally, for calculations we need the relation 

t a t a = C f l, (3.39) 

where Cf is the Casimir invariant of the SU(3) so that Tr t a t a = 4 in the funda- 
mental representation. 

The four quark operators we are dealing with are [9 

(qicr^qjq-jcr^qi) and (qa*P<h ]T q k jH a q k ) , (3.40) 

k=u,d,s 

for which the vacuum saturation gives 

4 4 
- ^(qiqt)(qjqj) and - -(<Mi) 2 , (3.41) 

so that finally the complete expressions read: 



D^f = (nii^mj) 2 (T^rriiigqiCrFqi) + m j (gq j aFq j ) ) (3.42) 
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where a is the coupling at the scale M 2 and (gqjaFqj) is the mixed condensate, 
for which we take (gcjiCrFqi) = 0.8(9^}. 

There is one last effect that we have to include in the OPE, namely instantons. 
M altman and Kambor have argued in Refs. [T3l l4*T] that these effects are impor- 
tant in (pseudo) scalar currents and should preferably be accounted for. This was 
done by taking a fit to several FESR and comparing the obtained result to the one 
from a BSR; it was found that consistency of the results strongly favors inclusion 
of instantons. Unfortunately, it is not entirely clear how to parametrize them and 
the model conventionally used in sum rules is the instanton liquid model (ILM), 
that goes back to Shuryak et. al. |42j . The OPE expressions for the FESR can 
be evaluated from 



\ f dss k U^(s) = T3K±r%) 2 %- J dss ^j lipiV - s)Yl ( Pl ^-s) , 



(3.43) 
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(3.44) 

where r] u d = 1, while r] us ps 0.6; Yi(s) and Ji{s) are the respective Bessel functions 
and pi is the average instanton size, which we will take to be pi ~ 1/0.6 as taken 
by Maltman and Kambor. In their calculations they add an additional error to 
their results, which they estimate from consistency. To us it is less obvious how 
to estimate this error, and we will not give one, especially since this uncertainty 
would not appear in our final result anyway. 



3.3 OPE for the Vector Contribution 



As we will see later, the OPE expressions we actually need for our calculations 
are the vector contributions. Leaving out dimension zero, the expression for 
dimension two is: 



D ij ( s ) 



(m} + mf) J2 c n (q 2 /p 2 )a n 



n=0 



Att 2 s 

±m, i m j e n (q 2 /p 2 )a n 



n=l 



n=2 



(3.45) 



As in the longitudinal contribution the quantity D L+T (s) will have to obey a 
RGE and the scaling of the coefficients will be according to the relations ()3.22j) . 
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The Qo, eio and fa are given in the appendix along with the d^. We are now 
in the position to give precisely the vector structure of the term 5^ from the 
previous section. Defining 
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(3.46) 



we get an expansion analogous to Eq. ()3.45|) : 
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(3.47) 



In writing down this line we have already added vector and axialvector terms 
$ij = ($ij,v + $ij a) which gets rid of the m^m.,- terms; since we discuss only the 
flavor symmetry breaking SRf the term proportional to /„(£) drops out, as do 
the contributions of m u , so that we finally arrive at 



5R 



M,L+T 



18S 



EW 



m 2 s (M 2 T 
Ml 



L+T(n) 
kl 



(3.48) 



n=0 



In accordance with Ref. ^U] we will refer to the sum in this equation as Ajy + . 

The dimension-4 structure of the vector current is more complicated than 
that of the scalar and we have therefore relegated a thorough discussion into the 
appendix and quote here the results. Adding vector and axial vector channels, 
the SU(3) breaking piece is given by 
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which leads to the following equation for the theoretical side of the r-decay rate: 
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(3.50) 



where all the coefficients and the definition of the expressions H kl ,kjt m 



L+T G L k + T and 

Qkf 1 are given in the appendix and we have defined the quark condensate oper- 
ator as 50 i = m s {ss — Eddd) . We emphasize that we have only made use of SU(3) 
symmetry and cancellation between axial and vector contributions and have not 
neglected any small terms, even though terms quartic in masses are expected to 
be negligible. We have, however, neglected the scaling of the quark condensate 
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operator as discussed before. For this operator we have taken the numerical val- 
ues for the quark condensates and as given in the appendix and estimated it 
as in Ref. PU]. 

Our values for H^ +T and G\f T as well as the Q^ +T are given in Table 13.21 
for all relevant moments. In principle one could use a combined fit for different 
moments to determine not only the strange mass, but also a more precise value 
of the quark condensates (for a thorough discussion of subtraction constants and 
numerical corrections to the GMOR relation see, for example, Ref. This is 

limited by experimental accuracy nowadays. 

The complete sum rule up to dimension-four is: 




Before we go on to combine the scalar and vector contributions let us discuss 
briefly the implications of dimension 6 contributions. We will have to include 
them in our analysis to be consistent since the leading dimension 6 terms turn out 
to be larger than the quartic mass terms of dimension 4. We expect the leading 
terms to be four quark operators that are not suppressed by quark masses and 
will estimate these along the lines of Ref. ^0] : if we take the OPE to be 

s 3 (D L J Tfi (s) - D L u r*(s)) = -350^^) £ q^fy?) a n (3-52) 

n=0 

and neglect the scaling of 4 quark condensate operators as for the regular quark 
condensates 3 we find for 5R kl 

5R k0 > 6 = c&R 00 ' 6 , (3.53) 

where the Cg are given by 

|z|=i 



3 The reason is somewhat different however: when discussing the renormalization group 
equations we have seen that the scaling is given by higher orders in the coupling as well as 
higher order logarithmic terms, both of which we do not include, so we must neglect the scaling 
altogether. The same argument holds for the scaling of the gluon condensate. 
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Then the final estimate for these terms is: 



«? 00 ' 6 - -^3 60^, (3.55) 
5R™>° ~ -^3 60^, 

1?7T 2 

rp20,6 2£Q L+T 
r 

Now we need to estimate the dimension 6 operators: again we will reduce them 
to the chiral condensates using the vacuum saturation hypothesis 



( , 36vr 2 +T 256tt 4 (ss) 2 - (dd) 2 
6R ~ -b EW ——dU 6 ~b EW ——a — 

T T 

~ 3SewC 6 ^^ « (0-6 ± 2.3) x lO" 3 , (3.56) 

where the numerical coefficients for the OPE are given in Ref. [HH] and we have 
taken the measured value of £°°' 6 = 0.001 ± 0.004 as in Ref. [TQJ and v s is the 
ratio of strange to non strange condensates, for which the current estimate is 

"■=£IH- 8±0 - 2 - < 3 - 57 > 

Our numerical result confirms the necessity to include these terms since the quar- 
tic mass terms are of order 10~ 5 and thus much smaller, as noted above. 

With all the tools given in this section it is generally possible to calculate the 
strange mass from the existing r-decay data. These determinations have been 
performed rather successfully "TT)[ ITT) l4~"""( 14*5] 147)] . but they all suffer from the 
same problem; to see this we shall very briefly go back to the perturbative series 
of dimension 2 and examine its convergence: we define 

Afc? = £ dn(Z) B u n) (3-58) 

n=0 

and 

Ai? = i(3Air ) + Ai?). (3-59) 
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Above we have stated that a contour improved calculation leads to better conver- 
gence and given the series for A^ 2 ) as an example. Now let us discuss the separate 
contributions. For the vector contribution we find 

A^ +T) = 0.7532 + 0.2137 + 0.0650 - 0.0002c 3 + ... (3.60) 
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(M) 

(0,0) 
(1,0) 
(2,0) 
(3,0) 
(4,0) 



A 



jl+tt 

kl 



0.98 ±0.13 ±0.10 
1.49 ±0.12 ±0.07 
2.02 ±0.25 ±0.03 
2.60 ±0.49 ±0.11 
3.23 ±0.78 ±0.21 



AW 
kl 

5.71 ±1.52 ±0.65 
5.73 ± 1.62 ±0.83 
3.98 ± 1.76 ±0.54 
2.48 ±2.05 ±0.65 
2.47 ± 1.59 ±0.35 



A 



2.16 ±1.22 ±0.09 
2.55 ±1.55 ±0.16 
2.37 ± 1.89 ±0.14 
1.94 ±2.23 ±0.20 
2.13 ± 1.88 ±0.13 



Table 3.1: Numerical values for the dimension 2 contributions. The first 
error shows the uncertainties for fixed coupling, while the second shows the 
uncertainty induced by the coupling constant. 



Taking C3 ~ C2 (02/01) = 323 as a naive estimate gives a last term of -0.0501, so 
that this series is rather well behaved. Unfortunately the same is not true for 
the longitudinal contribution. Even the improved series is unsatisfactory in its 
convergence: 

Aj? = 1.6334 + 1-2813 + 1-2892 + 1-5026 + ... (3.61) 

This bad convergence does not manifest itself too greatly in the central value 
of 5Rf because it is suppressed by a factor three when both contributions are 
added and the bad behavior is dominant only in higher orders. It does, however, 
have a significant impact on the errors obtained. Table l3~T1 lists our values for the 
dimension 2 contributions for the moments we will study. Throughout this work 
we will be using for a s (M T ) the value obtained by the ALEPH collaboration in 
Ref jTHj, meaning 

a s (M T ) = 0.334 ± 0.022 . (3.62) 

The origin of the errors is as follows: The theoretical error is taken by varying the 
parameter £ in an interval from 0.75 - 2.0. Additionally we have also included an 
error due to the uncertainty in the coupling constant. 

As mentioned before, the combined series is dominated in its convergence by 
the vector part of the sum for low orders. For higher orders the bad convergence 
of the longitudinal series causes the coefficients to increase again so that the series 
is merely asymptotic in some sense and we have to find a convention to truncate 
it. We have have followed the standard procedure to truncate the series after 
the minimum value coefficient and taken the full term of the next order as an 
additional error. To keep this discussion more transparent we will show exactly 
how the values for A kl arise and give the series explicitly 

A§j = 0.9732 + 0.4806 + 0.3710 + 0.3377 ± ... (3.63) 
Af } = 1.0396 + 0.5578 + 0.4829 + 0.4732 ± ... (3.64) 
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kl 



riL+'l ' 
U kl 
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kl 



(k,l) 



H, 



Q 



(0,0) -1.83 ±0.40 ±0.08 

(1,0) -1.34 ±0.75 ±0.20 

(2,0) -0.38 ± 1.01 ±0.30 

(3,0) 1.09 ± 1.10 ±0.39 

(4,0) 3.11 ± 1.36 ±0.44 



-1.30 ± 1.32 ±0.14 
-1.03 ± 1.76 ±0.28 
-0.43 ±2.02 ±0.42 
0.54 ± 1.96 ±0.53 
1.91 ±1.43 ±0.57 



0.076 ±0.047 ±0.162 
0.514 ±0.038 ±0.014 
0.932 ±0.030 ±0.008 
1.332 ±0.026 ±0.007 
1.714 ±0.077 ±0.024 



Table 3.2: Numerical values for the dimension 4 contributions; again we have 
split the uncertainty into the theoretical uncertainty and the uncertainty in the 
coupling. 



A$ = 1.1167 ±0.6443 ±0.6111 ±0.6382 + ... (3.65) 
A$ = 1.2008 ±0. 7393 ± 0.7569 ± 0.8361 ± ... (3.66) 
A$ = 1.2907 ± 0.8430 ± 0.9223 ± 1.0713 + ... (3.67) 

so that we will take the (0,0) and (1,0) series up to third order, the (2,0) up to 
second and truncate the (3,0) and (4,0) series after first order. For the longitudi- 
nal series itself one could also adopt this prescription but the behavior becomes 
extremely bad for higher moments, so that this does not seem meaningful. We 
have instead used an ad hoc definition and truncated the series at the same order 
as the total series. On the other hand we have seen that the vector series con- 
verges rather nicely and we will use the complete series with the estimated third 
order coefficient whenever we discuss Alf +T \ 
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Chapter 4 

Scalar and Pseudoscalar Spectral 
Functions 

You can know the name of the bird in all the languages 
of the world, but when you're finished, you'll know 
absolutely nothing whatsoever about the bird. 
(Richard P. Feynman, What Do You Care What 
Other People Think?) 

We have seen in the last chapter that and how the strange mass can be determined 
from r-decays. We have also seen that large uncertainties in this determination 
come from the bad convergence of the dimension two longitudinal perturbative 
series. This could be avoided by measuring longitudinal and transversal decays 
separately and then considering only the L + T terms of the sum rule. Un- 
fortunately, a complete spin disentanglement has not been performed over the 
entire allowed kinematic range. If the scalar and pseudoscalar spectral functions 
are known well enough up to an energy of M^, however, they can simply be 
subtracted from the decay rate using Eqs. ()3.8j) and (jrS.lOj) : 

«4 = -2ta» / i| (l - ^) (jj^J [/&'(.) -<£»(«)], I") 

where p^J(s) and p^J(s) are the longitudinal spectral functions for ud and us 
currents respectively. From Eq. ()3.6|) we see that 

pf{s) = (m + m,) 2 + (m* - m,) 2 , (4.2) 

which allows us to relate the scalar(s) and pseudoscalar (ps) spectral functions to 
the longitudinal ones we need (we will in fact relate them to those of the diver- 
gence of axial and vector currents, but this is equivalent, as shown in Eq. (j3.18|l ). 
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4.1 The Spectral Function for the ud Current 



As a general rule the ud current gives much smaller contributions to 5R than the 
us current and we could neglect the scalar channel entirely. This can be seen 
if we approximate the spectral functions by their theoretical expressions. The 
OPE of the pseudoscalar and scalar currents differ only in a general minus sign 
in the prefactor as well as some higher order terms that do not contribute much. 
The pseudoscalar spectral function is multiplied by (m^ + rrij) 2 , while the scalar 
comes with a factor of (m^ — rrij) 2 . This means the scalar spectral function for 
the ud current is suppressed by a factor of {m u — m^) 2 . Additionally, it is not 
entirely clear how to parametrize the scalar spectral function cleanly, so we take 
the theoretical expression as an estimate of how much the contribution might be. 

There is one more noticeable feature: dimensional analysis tells us that both 
spectral functions will rise linearly for large s. A resonance approximation can 
then only be justified for small s, where the theoretical expression is not yet a 
good approximation. This is to be compared to the vector sum rule that has a 
constant spectral density for large s, thereby making resonance saturation more 
intuitive. 

In any case the pseudoscalar channel gives sizable contribution and we will 
have to include it in our expressions. It is dominated by the pion pole and we 
will also add the next two higher lying resonances, namely the 7r(1300) and the 
7r(1800). In the narrow width approximation the spectral function is simply a 
sum of S functions, 

pi?(s)=E 2 /n^(s-M n 2 ), (4.3) 

where from the definition of the spectral function we see that 

(vr„,|9%|0) = v / 2/„M n . (4.4) 

Of course narrow width approximation is only good for the pion and we will 
include the excited resonances as Breit-Wigner functions: 



M 4 f 2 M 4 f 2 M 4 



where -Bj(s) is 



R i r,M, 

Bi{s) = -- 

71 



— T > ( 4 - 6 ) 
M 2 ) 2 + T 2 M 2 ] 

and 1,2 refers to the 7r(1300) and 7r(1800) respectively. The pion decay constant 
is known rather precisely to be f n = 92.4 ± 0.3 MeV, while the couplings for the 
higher resonances have been calculated by Maltman and Kambor in a combined 
fit of pseudoscalar sum rules They find 

fx = 2.20 ± 0.47 MeV , < f 2 < 0.37 MeV, (4.7) 

where we have added the errors they give in quadrature. 
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4.2 The Spectral Function for the us Current 



Since the contributions in the us channel are generally larger than those from 
the ud current, we have to include both scalar and pseudoscalar channels. The 
discussion of the pseudoscalar channel can proceed along the lines of the pseu- 
doscalar channel of the ud current. The three lowest lying resonances are the 
K meson itself, the if(1430) and the if (1830). As for the it a narrow width 
ansatz is again justified for the K, while for the excited resonances we need to 
use Breit-Wigner forms. The decay constants are 

f K = 0.113 ±0.001 MeV, f K (um) = 21.4 ± 2.8 MeV , (4.8) 



< /jr ( i830) < 8.9 MeV , 

that have been taken again from Maltman and Kambor 13: and the PDG. Nu- 
merical values for all resonance masses and widths are compiled in the appendix. 
There are rather large uncertainties attached to the parameters of the excited 
K resonances (in fact, both are listed by the PDG as "needing confirmation"), 
but these large uncertainties do not manifest themselves dramatically in the fi- 
nal values for SR™: varying the width of the If (1460) where the induced error 
should be most obvious gives an uncertainty of 0.000056, which is about 1% of 
the uncertainties in the coupling constants that give the dominant uncertainties 
as we will show in the next section. 

The scalar us channel could in principle also be described by a resonance 
ansatz as above, being dominated by the Kq(1430). The authors of jHj make such 
an approach by adding this as well as the next higher resonance, the _K"q(1930), 
and normalizing to threshold. Additionally the scalar form factor is included in 
an Omnes representation. The problems of this are twofold: on a simple level, 
trouble already arises when parametrizing the resonance (of course the same com- 
bined fit that was used for the pseudoscalar could be and has been used jl?] , but 
the uncertainties will be larger because there is no well known pole). Additionally 
there are significant dynamical contributions that interfere destructively in this 
channel and have to be included in a description. We will then have to think of 
something better, especially because the scalar us contributes significantly when 
adding up all four channels 

Following the approach for a dynamical ansatz by Refs. [JE1 EH], we take 
into account the lowest lying scalar states, which are excited Kit, Krj and Krj' 
resonances, and the phenomenological spectral function is given by 

3A 2 

P( s ) = ^f °k,{s)\F k ^s)\ 2 + a Kv {s)\F Kv {s)\ 2 + a Klll {s)\F Klll {s)\ 2 , (4.9) 
where A Kw = M\ - Ml, the okp{s) are the threshold phase space factors for a 
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Figure 4.1: The phenomenological spectral functions for the divergence of a 
vector current. We will use the central curve for our central values and take 
the others as error estimates. 



two body final state 



a KP (s) = 6 (s - (M K + M P ) 2 ) ^ 
and the form factors are defined by 



' x _ (M K + M P y \ / _ (M K - Mp) 



(4.10) 



(0\d,(s^u)\KP) 



' ^A Kn F KP (s) 



(4.11) 



These form factors have been calculated in Ref. [49J from KP scattering data, 
which allows us to give a phenomenological spectral function that should be valid 
up to about 2 GeV. It is important to note that all the spectral functions given 
above have already been written as spectral function for the r-decay current, 
while the scalar spectral function as plotted in Fig. 14.11 is specifically given for 
the divergence of the vector current. We will divide by the necessary factor of s 2 
in later calculations. We can now go on to subtract the pure longitudinal con- 
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tribution by evaluating the integral of (j4.1j) with the spectral functions discussed 
above. We shall begin with the pole contribution that can be directly seen to be 



5R L k f e 




2+k 



(4.12) 



where we have already taken / = 0. The integral for the remaining resonance 
contribution should only be taken from a threshold value s t hr to M 2 : 



-24vr z 



' as 



Ml 



s thr 



\ 2+k i v 

Ml) \Mf) 



(4.13) 



The threshold energy for the scalar us channel is Sthr = (Mk + M n ) 2 while for 
the pseudoscalar it is s t hr = (Mk + IM^) 2 . Equivalently, the threshold energy 
for the pseudoscalar ud current is 3M n . 



4.3 Pure (Pseudo) scalar Sum Rules 

In this section we shall study in detail the scalar /pseudoscalar part of the r- 
decay sum rule. To do so we will discuss separately the scalar and the pseu- 
doscalar channel of the ud and us currents and compare the OPE side with the 
phenomenological one. The results are given in Tables I4.1H4.5I They show the 
contributions SR^i (note the minus sign added in the strange channel) for the 
OPE as in 

R M = -4ni j ^Ti l {s/Ml)D^ L \s), (4.14) 

\s\=M* 

with the OPE discussed in the last chapter 1 and the phenomenological part as 
in Eqs. (I4.12j) and ()4.13j) . We have given all central values plus their respective 
errors and have additionally shown how the different terms in the OPE contribute. 
Clearly higher dimensional operators are strongly suppressed. 

We have used the following numerical values: m u = m s e u , m d = m s e d , 
(uu) = (dd) = -UM^/(m u + m A ) ■ (1 - 0.047), (ss)/(dd) = 0.8 ± 0.2 and 
(aFF) = (0.5 - 2) 2.1 • 10~ 2 GeV 4 . These are also the values that we will use in 
all subsequent calculations. Additionally we have taken m s (2GeV) = 105 MeV 
in accordance with the current world average to be able to provide a numerical 
value. The estimates are performed with the uncertainties given above, adding 

1 We find it more appropriate to truncate the perturbation series after the second order for 
all moments, which leads to better agreement between the expressions. 
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Phenomenology 


(1,0) 


-o.oo t to ± u.uuuuo ± u.uuuuo 






-U.UUlOZ IC U.UUUUO It u.uuuuo 


Pseudoscalar us 


(o,o) 


0.144 ± 0.024 trunc ± 0.009 scond 


Theory 


(1,0) 


0.140 ± 0.026 ± 0.009 




(2,0) 


0.138 ± 0.029 ± 0.009 


Scalar us 


(o,o) 


0.02790 ±0.0187 trMnc ± 0.0087 scond 


Theory 


(1,0) 


0.02758 ± 0.0200 ± 0.0088 




(2,0) 


r\ noon'? i r\ r\c\c\c\ i r\ r\r\oo 

0.02897 ± 0.0222 ± 0.0088 


Pseudoscalar ud 


(o,o) 


-0.0078 ±0.00015t rnnc ± 0.00002 mu ± 0.00005 md 


Theory 


(1,0) 


-0.0077 ± 0.00017 ± 0.00001 ± 0.00004 




(2,0) 


-0.0077 ± 0.00018 ± 0.00001 ± 0.00004 


Scalar ud 


(o,o) 


-1.74 ±1.28 trunc ± 0.14 md ± 0.32 mu • 10~ 5 


Theory 


(1,0) 


-1.73 ± 1.38 ± 0.14 ± 0.32 -HT 5 




(2,0) 


-1.81 ± 1.51 ± 0.14 ± 0.33 -lO -5 



Table 4.1: Total values and errors of 5R kl for all relevant currents and the 
three lowest moments; The errors are ordered as follows: For the phenomenol- 
ogy the dominant errors are from the coupling constants and we have given 
them in the order of the resonance, i.e. the first is for the it the next for the 
7r(1300) and the third (if given) for the 7r(1800). The scalar spectral function 
is varied in the limits given above. The theoretical uncertainty labled by trunc 
always includes scaling and truncation. 



the uncertainty that arrises by varying the parameter £ introduced in the last 
chapter from 0.75-2. We have also followed the standard procedure for asymp- 
totic series and taken the full value of the last coefficient as an error estimate, 
where we have truncated the series after the minimal coefficient. We have not 
taken into account the error induced by the strange quark mass for obvious rea- 
sons: we want to show that the errors for the theoretical expression are larger than 
for a phenomenological one and then result in larger uncertainties for a strange 
quark mass determination. It would thus be dishonest to include the error for 
the quantity we want to determine and pretend that this error plays a significant 
role. In any case the dominant uncertainties are those coming from truncation, 
scale dependence and the uncertainties in the quark condensates for the strange 
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(0,0) 


(1,0) 


(2,0) 


Dimension 2 


0.0503393 


0.0480009 


0.00476702 


Subtraction constant 


0.0829935 


0.0829935 


0.0829935 


Dimension 4 (mass) 


0.000738 


0.000983 


0.001247 


Dimension 6 (gluon condensate) 


-0.000395 


-0.000687 


-0.001024 


Dimension 6 (quark condensate) 


-0.000262 


-0.000480 


-0.000744 


Dimension 6 (mass) 


MO -8 


3-10 -8 


6-10" 8 


Dimension 8 


o i n — 5 

-3-10 


-0.000121 


-0.000277 


Instantons 


0.00885 


0.00770 


0.00603 


Total theory 


0.1436 


0.1395 


0.1376 


Phenomenology 


0.1345 


0.1210 


0.1102 



Table 4.2: Composition of the theoretical expressions for the pseudoscalar us 
current. 





(0,0) 


(1,0) 


(2,0) 


Dimension 2 


0.044778 


0.0426978 


0.0424037 


Subtraction constant 


-0.00869721 


-0.00869721 


-0.00869721 


Dimension 4 (mass) 


0.0007259 


0.0009666 


0.0012272 


Dimension 6 (gluon condensate) 


-0.000351 


-0.000611 


-0.000911 


Dimension 6 (quark condensate) 


0.000522 


0.000950 


0.001468 


Dimension 6 (mass) 


-MO- 8 


-3-10- 8 


-5-10- 8 


Dimension 8 


4-10" 5 


0.000144 


0.000330 


Instantons 


-0.009116 


-0.007871 


-0.00685 


Total theory 


0.02790 


0.02758 


0.02897 


Phenomenology 


0.02809 


0.01802 


0.01231 



Table 4.3: Composition of the theoretical expressions for the scalar us current. 



currents plus of course those from the quark masses in the ud currents. The 
phenomenological integrals are dominated so strongly by the poles (See Fig. 14 .2j) 
that the main uncertainties come from the tt and K coupling constants given 
by f w = 92 A ± 0.3MeV and f K = 113 ± O.OOlMeV. The uncertainties from the 
couplings to higher resonances are also significant, while those from masses and 
widths are negligible. 

In the tables given above there are two things that are especially noteworthy: 
first, that the central values for all currents agree rather well (though this agree- 
ment becomes worse for higher moments) and second, that the uncertainties of 
the phenomenological spectral functions are smaller than those of the theoretical 
one by an order of magnitude. This gives us much confidence in the general po- 
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(0,0) 


(1,0) 


(2,0) 


Dimension 2 


-0.000321 


-000306 


-0.000306 


Subtraction constant 


-0.007362 


-0.007362 


-0.007362 


Dimension 4 (mass) 


-5-10~ 9 


-7-10" 9 


-7-10^ 9 


Dimension 6 (gluon condensate) 


3-10~ 6 


4-10~ 6 


4-10~ 6 


Dimension 6 (quark condensate) 


1 1 c\ — 7 

1-10 


2-10 ' 


2-10 -7 


Dimension 6 (mass) 


-3-10- 14 


-7-10- 14 


-7-10- 14 


Dimension 8 


o i n— 7 
2-10 


6-10 


6-10 


Instantons 


-0.00011 


-9-10~ 5 


-8-10" 5 


Total theory 
Phenomenology 


-0.007789 
-0.007774 


-0.007757 
-0.007686 


-0.007739 
-0.007615 



Table 4.4: Composition of the theoretical expressions for the pseudoscalar ud 
current. 





(0,0) 


(1,0) 


(2,0) 


Dimension 2 


-2.658-10" 5 


-2.534-10" 5 


-2.517-10" 5 


Subtraction constant 











Dimension 4 (mass) 


-3-10~ 9 


-3-10- 9 


-4-10- 9 


Dimension 6 (gluon condensate) 


25 • 10~ 7 


4-10- 7 


5-10- 7 


Dimension 6 (quark condensate) 


-3-10- 8 


-5-10- 8 


-7-10" 8 


Dimension 6 (mass) 


2-HT 16 


5-10' 16 


9-10~ 16 


Dimension 8 


-2-10' 8 


-7-10- 8 


-2-10" 7 


Instantons 


9-10~ 6 


8-10- 6 


7-10- 6 


Total theory 


-1.740-10- 5 


-1.732-10- 5 


-1.810-10" 5 



Table 4.5: Composition of the theoretical expressions for the scalar ud current. 



tential of our ansatz. There is one more rather important fact that can be seen 
from those tables: a naive OPE estimate would suggest that the scalar spectral 
function should be very similar to the pseudoscalar because the OPE only differ 
by the global factor (rrii ±rrij) 2 , which is very similar for pseudoscalar and scalar 
if one quark is much heavier than the other. This similarity is obviously not 
present in the numerical values for the phenomenological spectral functions. The 
reason is that subtraction constants differ significantly in both channels (on the 
phenomenological side this is of course due to the fact that the weight functions 
put extreme emphasis on the low energy part of the spectrum and the main con- 
tributions come, as mentioned above, from the poles). 

Another check of the spectral ansatz can be to attempt an m s determination 
from pure scalar sum rules, and we will discuss this for our spectral ansatz. Let 
us study the divergence of vector current and denote the correlation function as 
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s[GeV 2 ] 



Figure 4.2: The k — weight function for r-decays (solid); due to the factor 
1/s 2 low energies are weighted very strongly. As a comparison we have plotted 
the Borel weight function (dashed) with a Borel parameter of M 2 =2 GeV. 



customary by 

^(q 2 ) = i I d'xe^(Q\T(f(x)f(0))\Q) , (4.15) 



with 

j(x) = df J ,q i j fl q j (x) = i{rrii - m^m^x) . (4.16) 

Clearly this quantity is well suited for such a determination because of the global 
mass factor. We can now generally attempt both Borel sum rules and finite 
energy sum rules, but we will concentrate on the scalar FESR since a Borel sum 
rule with the same scalar ansatz was performed by the authors of Ref. jT^ and 
pseudoscalar sum rules are discussed in much detail by Maltman and Kambor in 
Refs. for the ud and us currents. For finite energy sum rules it will be 

in interesting to vary the finite energy Sq. We can thus check different energy 
regions of the spectral function by using the stability of the obtained m s values 
as a criterion for the quality of our ansatz. We expect the results to be sensible 
only for a certain energy interval, since for low energies the OPE will break down, 
and for higher energies our spectral ansatz is now longer valid. 

If we now study finite energy sum rules for the correlator of axial vector diver- 
gences we see from dimensional analysis that we need two subtraction constants, 
i.e. two partial integrations, which leaves us with the FESR: 

1 S ° 
-— / dsW(s)V"(s) = [ dsw(s)p(s), (4.17) 
2m J J 

s=sq o 

where W(s) is the weight function integrated twice according to Eq. (|2.47jl . Since 
the pseudoscalar sum rules of Maltman and Kambor are rather well satisfied, we 
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will choose to work with their weight functions i.e. 



w(s) 



1- —) (l+A-)=>W(s 
s J V s J 



- (g ^\ a(s + s ) + 2s ), (4.18) 



12 



'0 



and set the parameter A = 0,2,4. Conveniently, Ref. |0] has already given some 
expressions for the second derivative and we need only quote the results. They 
are 



l + —a(-s) + 



+ 



/ 1995097 65869 



/5071 35 
715 



V 5184 



\ttu — mi 



216 
3 



+ 



1 40 / 
(rrii — mj) 2 



4tt 2 
miirij 



C(3)--C(4) + — C(5))a 3 (- S ) 
(m 2 + m 2 ) 



1 28 / 



C(3) a 2 (- S )(4.19) 



(4.20) 



— (aFF) + 

4 X 1 16tt 2 



+ 



4(m^ + rrij) — Sm^rn^- 

20 . ,n . . . .. 

ya(-s)j [rm{qiqi) +m j {q j q j )) 



+ 2 



25 



1 - ya(-s)J (miiqjqj} +m j (q i q i }) 



6(m; — mj 



-{mi(gqiaFqi) + m j (gq j aFq j }) 

4 4 4 

-47r 2 a(-s)-(g i g i )(g i g j ) - -vr 2 a(-s)-((g i g i ) 2 + (^gj) 2 ) 



(4.21) 



The expressions for dimension 4 and higher are not given explicitly in jS] , but can 
be easily calculated, while the O (a 3 ) coefficient was calculated from the compila- 
tion of scalar coefficients given in the appendix. Again, instantons are accounted 
for by using Eq. ()3.43|) . 



The strange quark mass determination will now proceed as follows: insert the 
spectral ansatz on the right hand side of Eq. ()4.17|) and plug in the OPE on the 
left. The resulting expressions can then be solved for m s for different values of 
so and the results are shown in Fig. 14.31 We see that the results are reasonable 
but only moderately stable for large A. The result for A=0 appears to be stable 
enough to possibly allow a determination, but this would obtain large uncertain- 
ties; we can show this by plotting the result for constant A and varying the scalar 
spectral functions in the limits given before. This is done in Figs. 14.41 and 14.51 
for A=0 and A=4. This results in a variation of m s of about ± 20, not including 
any uncertainties from the theoretical expression, a number which is clearly not 
competitive with other results. 
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0.2 



0.18 
0.16 



0.14 




0.08 



0.06 1 1 1 1 1 1 1 1 1 1 1 1 1 

2 2.5 3 3.5 4 4.5 5 

s [GeV 2 ] 

Figure 4.3: m 8 determination from pure scalar sum rules: we have plotted 
m s (2GcV) for A = O(solid), 2(dashed), 4(dot-dashed). 



Additionally we have performed the same calculations with the pseudoscalar 
sum rules discussed by Maltman and Kambor and the result is given in Fig. 14.61 
for comparison. Note that we have plotted the mass in the s region of 2-5 GeV 2 , 
while Maltman and Kambor discuss their sum rules only for s = 3-4 GeV 2 . 

As a conclusion, we hope to have shown that a FESR for pure scalar sum rules 
would be generally feasible but is only sensible with more precise scalar spectral 
data. Additionally, these plots serve as another confirmation of our spectral 
functions. 
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0.18 
0.16 
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m s [GeV] 
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0.08 
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Figure 4.4: m s determination from pure scalar sum rules: we have plotted 
m s (2GeV) for the three scalar spectral function given in Fig. 14. II for A=0. 
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Figure 4.5: As in Fig. 14.41 but for A=4. 
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Figure 4.6: As in Fig. 14. 31 but for the pseudoscalar sum rules. 
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Chapter 5 

Numerical Values for m s 



But if you 've ever worked with computers, you under- 
stand the disease - the delight in being able to see how 
much you can do. 
(Richard P. Feynman, 
Surely You're Joking, Mr. Feynman) 



The starting point of our numerical analysis is the vector side of the sum rule 
presented in chapter El Taking all terms up to dimension six, the complete 
equation reads: 



5R L + T 



12S 



EW 



3m 2 s (M?) 
2 Ml 



1 



el 



A 



(L+T) 



kl 



+ 3 



Mi 



el) 



(5.1) 



x 



H. 



L+T 
kl 



^e^fiki 



4tt- 



,50 4 (M} 
Mi 



-Q 



L+T 
kl 



+ 4 5R 00 > 6 , 



where 5R^ T is obtained from the ALEPH data of 5Ru quoted in Ref. [TT] as 
discussed in the last chapter by subtracting our ansatz for the scalar channels. 
From the tables in the last chapter we easily calculate the total values for SR^ 
and we give SR^ T , 6R& and 5R kl in Table O Since there is ongoing discussion 
about the size of \V US \ due to fact that the value obtained by Leutwyler and Roos 
50 from K decays shows a 2 o deviation from untarity, we have separated the 
uncertainties into those from the measured data and those induced by the CKM 



Matrix parameters; this uncertainty is dominated by the error in \V US 
absolute value is smaller. The central values for m», that we find, are: 



since its 





= 0.1973 GeV 




= 0.1640 GeV 


m 2 f 


= 0.1366 GeV 


m 3 f 


= 0.1152 GeV 


mf 


= 0.0986 GeV 
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(M) 




{[),[)) 
(1,0) 
(2,0) 
(3,0) 
(4,0) 


U.zli ±U.llo±U.UDz U.lOi ± U.UUO U.o ( 4 ± U. llo ± U.UDZ 

0.265 ±0.065 ±0.042 0.133 ± 0.003 0.398 ± 0.065 ± 0.042 
0.282 ±0.044 ±0.031 0.117 ±0.002 0.399 ± 0.044 ± 0.031 
0.292 ±0.034 ±0.024 0.104 ±0.002 0.396 ± 0.034 ± 0.024 
0.302 ±0.028 ±0.020 0.093 ± 0.002 0.395 ± 0.028 ± 0.020 



Table 5.1: r-decay widths decomposed into vector and scalar channels for the 
PDG2000 value of \V US \. 



Here we have taken the CKM Matrix elements from the PDG2000 unitarity fit 
as |K, | = 0.2225 ± 0.0021 and \V ud \ = 0.97525 ± 0.00046. We can perform the 
same calculations for the Leutwyler-Roos fit of \V US \ = 0.2196 ± 0.0026. Since we 
do not have access to the ALEPH data we will have to resort to the following 
procedure: the difference between the two values is 1.38 times the PDG error of 
| V us | and we can subtract from SRm the error induced by \ V US \ scaled accordingly 
to arrive at the values for 5Rm with the Leutwyler-Roos fit. The uncertainties 
induced by \V US \ were treated in a similar way, while we have not changed the 
given experimental uncertainties, since we do not know how much of it is actually 
due to the strange channel. 

The results for 5Rki and 5R^ T are given in Table IB~2*l the longitudinal con- 
tributions are obviously equal, so we can omit them. Again we calculate m s and 
arrive at significantly lower values: 





= 0.1508 GeV 




= 0.1415 GeV 




= 0.1225 GeV 


m 3,0 


= 0.1053 GeV 


mt'° 


= 0.0908 GeV 



It is now important to note the effects of replacing the scalar channels by a 
phenomenological ansatz in comparison to the OPE expressions: while the un- 
certainty in the longitudinal contribution could be suppressed to a point of neg- 
ligibility (see below), we have the drawback of being much more sensitive to 
experimental accuracy. This can be seen easily from the tables above and is par- 
ticularly drastic for the (0,0) moment: after we have subtracted the longitudinal 
contribution the experimental uncertainty is about one half of the absolute value 
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(M) 


SR^ 1 SRm 


(U, UJ 
(1,0) 
(2,0) 
(3,0) 
(4,0) 


(1 1 Q/( 4- n 1 1 9 4- fl D77 D OfiS 4- D 1 1 fi 4- D 077 
U.io4 ± U.llo ± U.lw I U.zoo ± U.llo ± U.Uf ( 

0.209 ± 0.065 ± 0.052 0.340 ± 0.065 ± 0.052 
0.242 ± 0.044 ± 0.038 0.356 ± 0.044 ± 0.038 
0.261 ± 0.034 ± 0.030 0.363 ± 0.034 ± 0.030 
0.276 ± 0.028 ± 0.025 0.367 ± 0.028 ± 0.025 



Table 5.2: Total and vector contributions to the T-decay widths for the 
Leutwyler-Roos value of \V US \. 

of 5Rqq T . From this observation we conclude that, with the present experimental 
accuracy, it should be worthwhile to include higher moments. These suffer from 
less experimental uncertainties since they are dominated by the better known low 
energy resonances. In contast, when using the theoretical expressions one usually 
restricts oneself to lower moments, since for higher moments the convergence of 
the OPE is not as good. In any case it is instructive to show this sensitivity by 
giving the results that arise, if SR^i is one experimental standard deviation below 
the central value, explicitly for the unitarity fit: 





= 0.130 GeV 




= 0.139 GeV 


m 2 f 


= 0.122 GeV 


m 3 8 >° 


= 0.105 GeV 


mf 


= 0.091 GeV 



These numbers are also in much better agreement with the results of other de- 
terminations and show less dependence on k, thus suggesting systematic experi- 
mental bias. If we now add up all uncertainties, we find for the final results: 









Exp. 




\v us \ 




Theor. 








m o,o = 


197 


± 


68 


± 


33 


± 


33 


± 


3 


MeV 


m]' = 


164 


± 


25 


± 


16 


± 


14 


± 


2 


MeV 


m 2,0 _ 


137 


± 


15 


± 


10 


± 


10 


± 


1 


MeV 


m 3,0 = 


115 


± 


10 


± 


7 


± 


11 


± 


1 


MeV 


m 4,0 = 


99 


± 


8 


± 


6 


± 


15 


± 


1 


MeV, 



where the theoretical uncertainties consist mainly of those from truncation and 
scaling. All others are rather small and listed in the table of App. A. 4, where 
we have also compiled all numerical input parameters. Uncertainties arising from 
scaling and truncation of the dimension 4 contributions are not listed, but are all 
smaller than 1 MeV. In our final results will include all uncertainties. 



53 



The weighting of these results has to take into account the correlations be- 
tween the different moments and we will take a conservative route in this work 
that should not overestimate the uncertainties too much, since these correlations 
are very large. We give the central value as the uncorrelated weighted average: 

m s = — ^Wim s>i . (5.2) 

^ i 

The weighting factors Wi are given by Wi = 1/af , with <7j the total uncertainty 
for m l, ° and w is just the sum over the In case of uncorrelated uncertainties 
the final uncertainty would simply be given by 1/ \/w, but in our case this leads 
to an underestimate. It should be conservative, however, to consider the smallest 
of the (Tj to be our uncertainty. 

This leads to a final value for the strange quark mass of 

m s (M T 2 ) = 122 ± 17 MeV. (5.3) 

At the standard scale of 2 GeV this gives 



m s (2GeV) = 117±17MeV. 



Performing all the error estimates with the Leutwyler-Roos fit, we find the fol- 
lowing results: 









Exp. 




\v us \ 




Theor. 








m o,o 


= 151 


± 


103 


± 


57 


± 


23 


± 


3 


MeV 


m]' 


= 142 


± 


30 


± 


24 


± 


12 


± 


1 


MeV 


m 2 s '° 


= 123 


± 


16 


± 


14 


± 


10 


± 


1 


MeV 


m 3,0 


= 105 


± 


11 


± 


10 


± 


13 


± 


1 


MeV 


4 


= 91 


± 


9 


± 


8 


± 


15 


± 


1 


MeV, 



and the central value turns out to be 

m s (M 2 ) = 107 ± 18 MeV . (5.4) 
Again we evolve to a scale of 2 GeV and find 



m s (2GeV) = 103±17MeV. 
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Chapter 6 

Conclusion, Discussion and 
Outlook 

I have to understand the world, you see. 

(Richard P. Feynman, 

Surely You're Joking, Mr. Feynman) 

We have calculated the strange quark mass from sum rules for r-decays and 
avoided the bad convergence in the scalar and pseudoscalar contribution to this 
sum rule by making a complete phenomenological ansatz for the scalar and pseu- 
doscalar spectral functions. This allows us to give not just an upper bound on 
the value, but to actually calculate it. 

As a comparison we show the results of other determinations in Fig. 16.11 
Except for earlier results from r-decays we give only the latest results. It seems 
that the results of different determinations have now finally become consistent, 
in contrast to the situation several years ago, when the results of, for example, 
sum rule and lattice calculations differed significantly. Note that our results too 
are in good agreement with the others. 

The values from r-decays are from Pich/Prades fU] and Chen et. al. 
Both use the total inclusive r-decay sum rule, but Pich/Prades also give an upper 
bound by making an ansatz for the spectral function that consists only of the kaon 
and pion pole and using the positivity of the spectral function. Chen et. al. take 
into account higher moments and use the same k dependent truncation that we 
discussed in chapter 3, while Pich/Prades use the full perturbative series for all 
moments. Nevertheless both agree very well. It is clear that even our conservative 
estimate has already reduced the uncertainties of the r-decay determination and 
we will discuss possibilities of further improvements below. 

The most promising other sum rule approach uses only the scalar or pseu- 
doscalar channels by considering the correlator of vector or axial vector diver- 
gences very similar to our final discussion in chapter |3J As mentioned, one can 
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Figure 6.1: m s (2GeV) for different determination methods. We have always 
quoted only the latest results, considering them an update of the previous ones; 
further references are given in the main text. 



use Borel sum rules [HI or finite energy sum rules J3] • in both cases the 

main uncertainty lies in the parametrization of spectral function, since higher 
K resonances are still rather poorly known and the spectral ansatz of ^2] for 
the scalar channel also leaves significant uncertainties. The pseudoscalar channel 
seems to allow for more precise determinations at the moment since it is strongly 
dominated by the kaon pole. In any case, significant improvements can only be 
made with experimental progress on the phenomenological spectral functions. It 
might, however, be interesting to use r-like weight functions, replacing by so 
and varying this parameter since the r-weights are sensitive to low energies and 
should effectively suppress anything but the poles. 

Another sum rule approach is the r-like e + e~ sum rule proposed by Narison 
|52j . Note that this has been criticised [H3] due to possibly large isospin breaking 
corrections. In its latest version [HI] that claims to be less sensitive to these 
effects one discusses the quantity 



3 IK 



All 



ud 



2-ko 2 



■s 



E\Y 



ds 11- 



1 + 



2s\ s 



(6.1) 
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as well the SU(3)-breaking combination 



Ai<J> = R T \ — Rr,<j) , 



(6.2) 



or isopin breaking 




(6.3) 



where R T) \ is the isovector and R T $ the isoscalar decay width. Just like the 
T-decay sum rules these sum rules suffer from cancellations that cause rather 
small a priori experimental uncertainties to become large. 

Finally, there are of course lattice results for m s , and we have quoted sepa- 
rately quenched [Uj and unquenched [53] values, since unquenched determina- 
tions tend to give smaller values and additionally the systematic uncertainty is 
unclear. On the other hand the quenched values are given by taking an average 
over all determinations and adding an uncertainty of ± 10 as systematic quench- 
ing uncertainty. It shall be very interesting to follow the progress here since the 
last word is most certainly not yet spoken. 

Going back to r-decays, we should ask ourselves about what further possibilities 
they hold. We have seen that our strange quark mass is extremely sensitive to 
experimental results, and with more precise data, as is expected from the coming 
charm and B factories, our results will improve significantly. This is not only 
due to purely experimental uncertainties: with better experimental precision we 
can restrict ourselves to lower fc-moments that also have a smaller theoretical 
uncertainty. Additionally, this should eventually get rid of the /c-dependence of 
the results, if its origin is indeed experimental, and then make the weighting a 
bit more reliable. Of course it will then become necessary to investigate in more 
detail the smaller uncertainties that have become irrelevant in our discussion. 

There is one last point that needs to be clarified, and that is the question 
of the CKM element \V US \. We have shown that due to the cancellations that 
make the experimental uncertainties large, we are also rather sensitive to the 
value of this parameter. This can result in one of the following conclusions: first, 
that we can improve our strange mass determination if this question is answered 
and second that we can (as has been done before) turn our investigation around 
and provide input for m s in order to determine 5R and from that try to obtain 
values for \V US \. At the moment this is of course somewhat speculative because 
it relies very much on further experimental input, but the possibility should not 
be entirely neglected 1 . 



The potential of this approach is discussed in detail in Ref. [HI] 
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Appendix A 

Miscellaneous Coefficients and 
Constants 



A.l Coefficients for Running Masses and Cou- 
plings 

There are multiple conventions for defining the (3 function, that mostly differ 
by the sign and/or a factor 2, that arises if the derivative is not taken to be 
with respect to \i but to /i 2 . We have adopted the convention of Pascual and 
Tarrach [2B], except for a minus sign and an extra factor a, so that the expansion 
coefficients for the function in the MS scheme are 

0i = | (11JV - 2/) , /3 2 = ±(l7N 2 -5Nf-3C f f) , (A.l) 
1 /2857. 1415 „ , 79 2 205 ^ f , 11 f2 , ^ 2 



03 = — N a N z f + —Nf NCff + —C f f z + Cif 

Hi 32 V 54 54 J 54 J 18 tJ 9 tJ fJ 

140599 445 

/3 4 = + C(3) • 

^ 2304 16 w 

For 7(a) we use 

7i = \c f , ^=2^(97 # + 9(7,-10/), (A.2) 
C, film at2 129 N(Jf (278 , o ^ „ , , 129 ^ 2 

35 



73 = AT NCf-[ + 24C(3) JV/+ C f 

16 32 V 108 4 ' V 27 sv V 2 7 



23 - 24C(3)jC / /-^/ 2 

2977517 9295 _ 135.,,, 125 . /rN 
^ = ^073^-W« 3) + - C(4) --« 5) - 

Here A" is the number of colors, / the number of flavors and Cf is the Casimir 
invariant, Cf = (N 2 — 1)/2N. For simplicity we have given both $4 and 74 for 
three flavors and colors only. 
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A. 2 Weight Functions for Higher Moments 

A general representation of the Function \ where j is either L or L + T, is 



x) 



2{i-xf +k Y: 



/! 



(x -1) 



3(1 -x 



,3+fc 



i\ (x - iy 



(6 + k + n) + 2(3 + k + n)x 
(3 + k + n)(4 + k + n) 



E 



n =o _ n )' n ' ^ + k + n ' 
Additionally, we give the functions for the values of k and I that we consider: 



(k,l) 



(0,0) 
(1,0) 
(2,0) 
(3,0) 
(4,0) 



x 



{l-xY{l + x) 



_L(l_ x )4 (7 + 8x ) 
^(1 -x) 5 (4 + 5x) 



^(l-x) 6 (9 + 12a;) 
^(l-x) 7 (10 + 14x) 



X 



!(i 



1(1 -xf 



|(1 -.x 

id 



A. 3 Resonance Parameters 

For the pseudoscalar spectral functions we need the following parameters of low 
lying resonances, that have been taken from |T3J |2H E] • 





Mass(MeV) 


Width(MeV) 


Decay constant ( Me V) 


7T 


138 




92.4 ±0.3 


tt(1300) 


1300 ± 100 


400 ± 200 


2.20 ±0.46 


tt(1800) 


1800 ± 13 


210 ± 15 


0.19 ±0.19 


K 


495 




113 ±1 


K(U60) 


1460 ± 60 


260 ± 10 


21.4 ± 2.8 


AT(1830) 


1830 


250 


4.5 ±4.5 



The uncertainty of the 7r(1300) width and coupling are correlated, since some 
of the uncertainty from the width manifests itself in the determination of the 
coupling constants, but this correlation is not important for our strange mass 
determination. Since the excited K resonances are listed by the PDG as needing 
confirmation, they do not give uncertainties for widths and masses. They do, 
however, give several central values (in the case of the i^(1460)), and the numbers 
we quote are estimated from them. 
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A. 4 Numerical Input Parameters 



For reading convenience we compile here all our numerical input parameters in- 
cluding the corresponding uncertainties. Finally, we give the induced uncertain- 
ties for m s . 



Parameter 


Value 


AmJMeV) 


« s (M r 2 ) 


0.334 ±0.022 


2 


(Ss) 


0.8 ±0.2 (uu) 


7 


m u 


m s (0.029± 0.003) 


o(io- 3 ) 


m d 


m s (0.053± 0.002) 


^(io- 1 ) 


(uu) 


-/A/K + %)(1 - 0.047) 


see m u and m d 


Sew 


1.0201 ±0.0003 


< c(io- 3 ) 


(aFF) 


(0.5-2)2.110- 2 GeV 4 


< ^(io- 3 ) 1 


Vud 


1 






0.6 




Pi 


1/0.6 




(gqiaFqi) 


0.8 (q iqi ) 





According to the way we estimate the uncertainties, we give here those un- 
certainties that contribute to the moment with the smallest total uncertainty. 
Additionally we need the CKM matrix elements \V US \ and \V ud \. For \V ud \ we 
take 0.97525 ± 0.00046, while for \V US \ we discuss both the PDG unitarity fit 
\V US \ = 0.2225 ± 0.0021 and the Leutwyler/Roos value of \V US \ = 0.2196 ±0.0026. 
Uncertainties for the strange mass determination are discussed in the main part 
of the text. 



1 In our calculations the gluon condensate drops out, but it appears in higher dimension 
contributions suppressed by quark masses. 
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Appendix B 



Higher dimension terms for the 
OPE 



B.l Dimension-two Coefficients 

The needed coefficients for the longitudinal contribution are: 

17 , 9631 35 A/ , .„ , 

rf 00 = l, d 10 = —, d 20 = - — C(3) , (B.l) 

4748953 91519 A/oN 715 A/rN 

d 30 = C(3) C(4) H C(5) • 

30 5184 216 SV ; 2 SV ; 12 sv ; 

For the transverse: 

179 520 

coo = l, cio = ^, c 20 = + — C(3) - — C(5) , (B.2) 

5 

eoo = , eio = ^ , e 20 = ^ - ^C(3) - ^C(5) , 



/oo = 0, /; 



10 



13 




23077 




C20 = 




T ' 


432 + 


2 




769 55 


3' 


e2o = 


~54~ ~ 27 




/20 = 


32 8, 


0. 


-¥ + 3< 



The notation is mainly according to Ref. JO], albeit slightly different, where all 
of the coefficients have been compiled. 



B.2 Dimension- four Coefficients 

We will present here the complete L + T dimension 4 contribution. This includes 
operators m 4 , quark condensates multiplied by masses and finally the gluon con- 
densate operator, so that the general expression reads: 
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where 

^n +T (s/V) = \{aFF)p L n +T (s/^) + Y,{rn k Mk)r L n +T (s/ l ?) (B.4) 
6 T 

+2{m i q i q i + m j q j q j )q^ +T {s/ > 2 ) - — (mf + mfj h^ +T (s/fi 2 ) 
±g(™35i5i + m i ^^)^ +T (s//x 2 ) + ^m 2 m 2 g^ +T (s/ /i 2 ) 

(m 2 + m 2 ) ^(s//* 2 ) - JiE^W) 

k 



71 W 



This simplifies considerably when the only SU(3) breaking contributions are 
taken, and we have tried to write the equation in such a way as that these 
simplifications become most easy to see. We will refer to the terms by the letters 
of the coefficients that they carry. Obviously the first line will drop out, as will 
the u and the j terms in the last two lines. Additionally the t and k terms drop 
out if we add vector and axial vector channels. The remaining terms can be 
written as 

* 2 {Duf A (s) - D L u f>\s)) = ~A50 4 (-s) E rf/zi 2 ) + (B.5) 

n 

^rnt (l - e 2 d ) £ [(l + e 2 d ) h? T (s/f) - e 2 u g L n + T (s / » 2 )} a" . (B.6) 

" n 

We can now easily arrive at Eq. ()3.49|) . if we introduce the following notations: 
Qtt T = ^E?n +T (0 / ^ k L l +T (s/M> n (-es) , (B.7) 



\ 8 \=M* 



4 



n \s\=M* V \ T) / 



4 



s|=A/| 

We have set /i 2 /s = — £ 2 to allow for an error estimate. The parameter £ is later 
varied in the interval 0.75 < £ < 2. 
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The C n o have been calculated to second order for the condensates and, as given 
in [TU|, are 



P L *r = o, P L l0 +T = i, P L 2 r 



r-L+T _ n L+T ... n L+T 

00 — u ; '10 — u ) '20 



?ob +T — 1 ) Qio T — ~ 1 j <?2o +T 



_ n + L+T _ n -iL+T 

^00 — u ; fc 10 — 1 5 t 20 



7 

6' 

131 



24 ' 



17 

y 



We need scaling of these coefficients as an error estimate, but to the order we are 
calculating there are only terms linear in the logarithms. The nonzero ones are: 

r L+T _ _ - n L+T _ 9 + L+T _ 9 

'21 — g ; H21 _ 4' 21 ~~ 4' 

On the other hand the mass coefficients have only been calculated up to first 
order, which is, including the scaling behavior: 



h L o +T = h L S T = | - 2C(3) - |L + L 2 



k L , +T = 0, k^ T = 1-h, 

5 

g L * +T = i, sf +T = y-^C(3)-4L, 

j i+T = o, = o, 

4 +T = o, u L x +T = o. 

2 

where L = hi^-. The longitudinal coefficients are 

41 14 
Jtf = fcf = --2C(3)- y L + L 2 , 

= fcf = 8-^(3)-yL + ^L 2 

Jo" = 0, = 0. 
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Appendix C 



Proof of the Formula for 
r-Decays 

Know how to solve every problem that has ever been 
solved. 

(Richard P. Feynman) 

The central formula of this work (this is by no means to suggest that it is new, 
but all that we are doing relies on it) is ()3.8)1 that we have stated without proof. It 
was originally derived by Tsai in 1971 [HE] in the framework of Current Algebras. 
In this chapter we present a somewhat more elaborate and different approach 
that basically relies on the optical theorem and the Cutkosky cutting rules that 
follow from it. In any case of doubt we will be using the notation of Peskin and 
Schroeder [22] . 

First, remember that we have decomposed the correlator as follows: 

W v (q) = {<fq v - q 2 g^)lf r \q 2 ) + q^q u U (L \q 2 ) . (C.l) 

We shall first be concerned with the transversal part, and we will see that the lon- 
gitudinal component can then be easily inferred by the calculation. The general 
formula for a decay rate T is: 

r -2k! (n (g^) - W*o» - x» . (c.2) 

where \A4\ 2 is the squared decay matrix element, and we have let a particle A 
decay into final states /. 

According to the optical theorem this squared matrix element as in Fig. IC.ll b) 
is equal to twice the imaginary part of the doubled diagram of Fig. IC.ll a) - these 
are the aforementioned Cutkosky Rules. In this process a we have to sum over 
all possible hadronic final states and integrate over their respective phase space 
factors. Then the delta function disappears and only the phase space integration 
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a) b) 



Figure C.l: The squared matrix element for r decays into hadrons b) can be 
calculated as the imaginary part of the connected diagram a); a phase space 
integral as well as a sum over all possible hadronic states is understood. 



over the momentum of the neutrino is left. First, however, we need to calculate 
the matrix element for Fig. IC.ll a). which is: 

.2 'l-y^ 



iM 



-i) 2 |-w(p T )7 M 



2 J 

'1 -7 5 ' 



M 2 W 



I u(p T 



M 2 W 14 2 



(C.3) 



Here we have used the subscript v for a Lorentz index and to denote the momen- 
tum of the neutrino. It should be obvious from the context, however, which of 
the two is meant in all steps. 

Let us discuss the treatment of the spinors now and for that drop all the 
irrelevant constants. We need to sum over the final and average over the initial 
spin: 



2 12 U (Pr)l^ 
spins 



'1 -7 s 



2 



r7M 



'1 -7 5 



u(p u )u(p v )j, 
/ 1 -7 5 ' 



Ufa) 



1, 



-Tr 



r7M Pulu 



1 -7 5 



(C.4) 

(C.5) 
(C.6) 
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= 0« + - gTto ■ Vr - ie^PraPvp) • (C.7) 

This expression has to be dotted into the matrix element (|C.3J) . In doing so the 
£ afi/3u ^ erm w j}] d r0 p out due to symmetry and we obtain: 

^El^l 2 = i ^r(«+J'X-?^-p T )(w,-«v) n(r) (*) (C8) 

spins W 

[sp + -p T + 2(p + q)(qp T )}U {T) (s). 



16M& 

In the case of massless neutrinos 1 

2p T p u = 2p v q = M 2 T - s, 2p T q = M 2 T + s , (C.9) 

so that the after taking the imaginary part and multiplying with the necessary 
numerical factors our Matrix element is: 



(CIO) 



Finally we need to multiply with the phase space factor. Changing the integration 
variable from p v to s we see from 

that dp u = —ds/2M T . Additionally for the neutrino p = E and our total phase 
space factor turns out to be 

and the minus sign vanishes if we choose the appropriate order of integration 
bounds, i.e. integrate from to M^. Throwing all of this together and performing 
the angular integral the decay rate is 

M 2 2 



2M T 

At this point we can go back to include the longitudinal part. The calculation 
is in principle similar, except for Eq. (jC.8|) . If we take only the q^q v term the 
second line becomes 

sp + -p T + 2(p + q)(qp T )}rt L \s). (C.14) 



16M, 



1 The following relations are most easy to see in the CMS 
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Adding this to the result for the transverse correlator we arrive at the final result 
for the decay rate of the r into hadrons: 



r(r — > hadrons) 



1 7 ds g 4 ( s \ 

1 + 2j^J Imn (T) (s) + Imn (L) (s) 



The decay rate of the r into leptons is 



r(r -> ez/ r z/ e ) 



192,3 ' < ai6 > 
where the Fermi constant GV is defined as Gf/V% = g 2 /(SM^), so that we find 



r(r -> ez/ T z/ e ) 



<? 4 M r 5 



32 • 192M^tt 3 ' 



(C.17) 



Finally we can now give the expression for R T 



Ml 



R T = 127T 



ds ( s 



1 + 2— j imIT (T) (s) + Imn (L) (s) 



(C.18) 



which is the desired result. 
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